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the  user  to  specify  the  degradation  behavior  for  each  mode  of  failure. 
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SUMMARY 


Impulsive  loading  of  structures  composed  of  fiber-composite  materials  has  appli¬ 
cation  in  many  applied  fields.  Tests  on  fiber-composite  motor  case  models  indicate 
a  wide  range  of  material  damage,  including  matrix  cracking,  delamination,  and  fiber 
breakage  resulting  in  catastrophic  failure  of  the  case.  A  composite  damage  model 
was  recently  implemented  into  the  nonlinear  finite  element  code  DYNA3D  (Refer¬ 
ence  [1])  to  analyze  impulsively  loaded  materials  and  structures.  The  validity  of  the 
composite  damage  m.odel  is  assessed  in  the  present  report  by  comparing  the  model 
with  experimental  results. 

This  summary  and  the  main  body  of  the  report  are  in  two  sections.  The  first 
section  is  a  review  of  commonly  used  failure  criteria  and  the  development  of  the 
Chang  interactive  criterion  currently  implemented  in  DYNA3D.  The  second  section 
compares  failure  criteria  and  post-failure  degradation  rules  with  experimental  results 
to  assess  the  validity  of  the  DYNA3D  composite  damage  model.  Modifications  to 
both  the  failure  criteria  and  degradation  rules  are  recommended. 

Review  and  Development  of  Failure  Criteria. 

Failure  criteria  are  used  to  relate  critical  combinations  of  stress  or  strain  to  failure  in 
a  material.  Four  commonly  used  failure  criteria  are  the  maximum  strain,  maximum 
stress,  Tsai-Hill,  and  Tsai-Wu  criteria.  These  criteria  are  all  in  agreement  for  the 
uniaxial  failure  stresses  in  the  principal  material  directions  of  a  lamina.  The  criteria 
disagree  on  what  constitutes  failure  for  biaxial  stress  states. 

The  maximum  stress  and  maximum  strain  criteria  predict  five  independent  failure 
modes;  tensile  and  compressive  longitudinal  failure,  tensile  and  compressive  trans¬ 
verse  failure,  and  shear  failure.  The  Tsai-Hill  and  Tsai-Wu  criteria  predict  only  the 
onset  of  failure,  not  the  mode  of  failure,  and  are  called  interactive  criteria  because 
failure  depends  on  more  than  one  stress  component. 

Hashin  in  Reference  [2]  developed  a  set  of  interactive  failure  criteria  in  which  dis¬ 
tinct  fiber  and  matrix  failure  modes  are  modeled.  Chang  in  References  [3-5]  modified 
Hashin’s  criteria  to  include  nonlinear  shear  behavior.  Three  failure  modes  are  mod¬ 
eled:  fiber  breakage,  matrix  cracking,  and  matrix  compression  failure.  Fiber  compres¬ 
sive  failure  is  not  modeled.  The  Chang  criteria  are  implemented  in  DYNA3D.  When 
failure  occurs,  some  or  all  of  a  lamina’s  elcistic  constants  are  set  to  zero,  depending 
on  the  mode  of  failure. 

It  is  important  to  identify  the  fiber  and  matrix  failure  modes  because  the  effect 
of  each  mode  on  the  ultimate  strength  of  a  layered  composite  is  quite  different.  For 
example,  non-catastrophic  matrix  failure  generally  precedes  catastrophic  fiber  failure 
in  filament-wound  motor  cases.  Some  factors  which  affect  the  fiber  and  matrix  failure 
modes  are;  properties  of  the  fiber  and  matrix,  fiber  volume  fraction,  and  the  fiber- 
matrix  bond  strength. 
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Verification  of  DYNA3D  Composite  Damage  Model. 

Calculations  were  performed  to  verify  the  implementation  of  the  failure  criteria  and 
degradation  rules  used  in  the  DYNA3D  composite  damage  model.  Each  failure  mode 
and  degradation  rule  was  demonstrated  through  calculations  on  single-element  shell 
laminates.  We  are  satisfied  that  the  failure  criteria  and  degradations  rules  in  our 
DYNA3D  version  are  implemented  correctly. 

Computational  results  are  also  compared  with  experimental  results  found  in  the 
literature  to  establish  the  validity  of  the  failure  criteria  and  post-failure  degradation 
rules.  The  Chang  failure  criteria  are  compared  with  three  sets  of  material  property 
tests: 

•  lamina  off-ajcis  tests  to  verify  the  tensile  failure  criteria  (fiber  and  matrix). 

•  lamina  combined  stress  tests  to  verify  the  matrix  failure  criteria, 

•  biaxial  stress  tests  on  quasi-isotropic  laminates  to  verify  the  fiber  failure  criteria. 

The  Chang  fiber  breakage  criterion  is  in  agreement  with  biaxial  stress  test  data  if 
reasonable  assumptions  are  made  regarding  tnatrix  cracking  prior  to  fiber  breakage. 

However,  the  comparison’s  indicate  that  the  Chang  matrix  failure  criteria  do  not 
adequately  model  matrix  strength  under  combined  loading  conditions.  A  modified 
form  of  the  Tsai-Wu  criterion  is  in  better  agreement  with  the  test  data  and  is  therefore 
recommended  for  implementation  into  DYNA3D. 

Most  analytical  and  experimental  results  available  from  the  literature  determine 
the  post-failure  stress-strain  response  of  laminates  with  cracked  plies.  One  common 
method  for  estimating  stiffness  reductions  in  laminates  with  cracked  plies  is  the  ply 
discount  method.  With  this  method,  one  or  more  of  the  eleistic  constants  of  a  cracked 
ply  are  set  equal  to  zero.  A  ply  discount  scheme  is  currently  implemented  in  DYNA3D. 

Tests  and  analyses  results  indicate  that  the  reduction  in  ply  properties  due  to 
matrix  cracking  depends  on  the  ply  stacking  sequence,  ply  stiffnesses,  and  cracking 
ply  thickness.  The  ply  discount  scheme  does  not  adequately  predict  the  reduction  in 
stiffness.  Thus,  a  general  post-failure  degradation  rule  should  be  implemented  into 
DYNA3D  that  is  capable  of  modeling  ply  stiffness  reductions  for  arbitrarily  stacked 
laminates.  A  flexible  post-failure  degradation  rule  is  recommended  and  modifications  » 

to  the  current  implementation  are  suggested  which  allow  the  user  to  specify  the 
degradation  behavior  for  each  mode  of  failure. 
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PREFACE 


This  topical  report  is  part  of  an  effort  to  develop  and  verify  numerical  techniques 
for  modeling  fiber-composite  materials.  Numerical  results  are  compared  with  exper¬ 
imental  results  available  from  the  literature  to  test  the  validity  the  fiber-composite 
failure  criteria  and  post-failure  degradation  rules  implemented  in  DYNA3D.  Numer¬ 
ical  results  which  verify  the  constitutive  behavior  of  the  laminated  shell  model  were 
included  in  a  previous  report,  Fibrous  Composite  Laminates:  An  Introduction  to  the 
Theory,  and  Verification  of  the  DYNA3D  Implementation,  DNA-TR-88-191. 

The  DNA  technical  monitors  were  Dr.  Michael  Frankel  and  Maj.  Paul  J.  Wolf. 
The  APTEK  program  manager  was  Mr.  Stephen  H.  Sutherland.  Suggestions  and 
comments  on  the  topical  report  were  provided  by  Drs.  Leonard  E.  Schwer  and  Herbert 
E.  Lindberg  of  APTEK. 
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SECTION  1 

BRIEF  BACKGROUND  ON  THE 
DEVELOPMENT  OF  COMPOSITE  DAMAGE 

MODELS. 


Failure  criteria  relate  critical  combinations  of  stresses  or  strains  to  failure  in  a  mate¬ 
rial.  Some  criteria  predict  only  the  onset  of  failure,  while  more  sophisticated  criteria 
predict  the  mode  of  failure  as  well.  VVe  begin  this  section  by  reviewing  some  of  the 
failure  modes  observed  in  fiber-composite  materials.  The  subsection  on  failure  modes 
is  followed  by  a  review  of  failure  criteria  that  are  commonly  used  by  analysts  today. 
Finally,  the  development  of  the  Chang  composite  damage  model,  which  is  currently 
implemented  in  DYNA3D,  is  presented. 

1.1  FAILURE  MODES  IN  UNIDIRECTIONAL 
COMPOSITES. 

Distinct  fiber  and  matrix  failure  modes  are  observed  in  fiber-composites.  It  is  impor¬ 
tant  to  distinguish  the  modes  of  failure  because  the  effect  of  each  mode  on  the  ultimate 
strength  of  a  layered  composite  is  quite  different.  For  example,  non-catastrophic  ma¬ 
trix  cracking  generally  preceeds  fiber  breakage  (catastrophic)  in  filament-wound  motor 
cases.  Matrix  cracking  of  at-angle  plies  leads  to  load  redistribution  between  the  hoop 
and  at-angle  layers,  while  fiber  breakage  often  leads  to  catastrophic  failure  of  the  mo¬ 
tor  case.  The  factors  affecting  the  fiber  and  matrix  failure  modes  are  briefly  discussed 
in  this  subsection.  Most  of  this  information  was  obtained  from  the  lecture  notes  of  a 
UCLA  short  course  on  advanced  analysis  of  composite  materials,  Reference  [6].  The 
course  notes  on  failure  modes  were  presented  by  Mr.  Longin  Greszczuk  of  McDonnell 
Douglas  Space  Systems  Company  (MDSSC). 

1.1.1  Fiber  Tensile  Failure. 

Some  factors  which  affect  tensile  failure  in  the  longitudinal-fiber  direction  are: 

•  properties  of  the  fiber  and  matrix, 

•  fiber  volume  fraction, 

•  fiber-matrix  bond  strength, 

•  fiber  length, 

•  nonuniformities  (scatter)  in  fiber  strength,  modulus,  or  pre-tension. 
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Figure  1.  Tensile  failure  mode  in  longitudinal  fiber-direction. 


Fiber  misalignment  and  test  methods  can  also  affect  the  measured  strength.  For 
example,  a  graphite-epoxy  specimen  tested  with  fibers  oriented  at  4°  to  the  test  axis 
exhibits  a  drop  in  measured  strength  of  about  17%.  This  is  because  the  biaxial 
stress  field,  in  the  principal  material  directions  of  the  composite,  causes  fiber-matrix 
debonding  or  matrix  cracking  failure  at  the  low  transverse  tensile  and  shear  strengths 
of  the  material.  Matrix  failure  reduces  the  ability  of  the  composite  to  carry  additional 
load  in  the  longitudinal-fiber  direction.  The  ability  of  the  matrix  to  transfer  load 
by  shear  to  an  unbroken  section  of  fiber  is  demonstrated  by  comparing  the  ultimate 
strength  of  fiber  bundles  tested  with  and  without  a  matrix.  Greszczuk  in  Reference  [6] 
reports  that  the  ultimate  strength  of  fiber  bundles  tested  without  a  matrix  material 
is  up  to  40%  less  than  that  of  a  fiber-matrix  composite. 

Tensile  tests  in  the  fiber  direction  of  unidirectional  polymer- matrix  specimens 
indicate  that  the  failure  surface  is  jagged,  discussed  in  References  [7]  and  [8]. 
Random  fiber  breaks  lead  to  matrix  cracks  both  transverse  to  the  fibers  and  along 
the  fiber  direction,  as  shown  schematically  in  Figure  1.  Transverse  cracks  coalesce, 
reducing  the  load  carrying  capacity  of  the  specimen,  until  complete  fiber  breakage 
failure  occurs. 

Longitudinal  tensile  strength  is  commonly  modeled  with  either  a  rule  of  mixtures 
or  a  statistical  flaw  model.  The  rule  of  mixtures  relates  the  longitudinal  tensile 
strength  of  the  composite  to  the  strength  of  the  fiber  and  matrix  components  in 
proportion  to  the  volume  fraction  of  each  component.  The  rule  of  mixtures  tends  to 
overestimate  the  strength  of  composites  with  poor  fiber-matrix  interfaces.  The  flaw 
model  considers  the  statistical  nature  of  fiber  breakage  by  including  a  constant  which 
is  a  mecisure  of  scatter  in  fiber  strength.  Fiber  length  is  included  in  determining  the 
scatter  because  longer  fibers  have  more  flaws,  hence  break  at  lower  failure  stresses 


than  shorter  fibers.  An  ineffective  fiber  length  is  also  included  in  the  model.  When  a 
fiber  breaks,  the  load  is  transferred  by  shear  in  the  matrix  back  to  the  fiber  a  short 
distance  along  its  length  from  the  break.  The  ineffective  fiber  length  is  that  length 
required  to  pick  up  a  certain  percentage  of  the  load,  and  is  thus  a  measure  of  the  load 
redistribution  capability  of  the  composite.  Further  discussion  of  the  statistical  flaw 
model  is  given  in  References  [8-9]  by  Rosen. 

1.1.2  Fiber  Compressive  Failure. 

Compressive  failure  in  the  longitudinal-fiber  direction  is  strongly  influenced  by  the 
following  factors: 

•  fiber  and  matrix  properties, 

•  fiber-matrix  bond  strength, 

•  the  fiber  volume  fraction. 

Possible  failure  modes  include  microbuckling  of  bonded  or  debonded  fibers,  inter¬ 
action  failure,  or  direct  fiber  failure.  These  failure  modes  are  reported  by  Greszczuk 
in  Reference  [10]  and  discussed  in  the  following  paragraphs. 

Microbuckling  with  Bonded  Fibers.  For  fibers  well  bonded  to  the  matrix,  this 
type  of  microstructural  failure  can  occur  under  compression  in  either  the  extensional 
mode  or  shear  mode,  as  shown  in  Figure  2  and  discussed  in  References  [1 1]  and  [12].  In 
the  extensional  mode,  the  matrix  extends  or  contracts  in  the  direction  perpendicular 
to  the  fibers.  In  the  shear  mode,  the  matrix  deforms  primarily  in  a  shear  deformation 
between  adjacent  fibers.  The  mode  depends  on  the  fiber  volume  fraction,  and  moduli 
of  the  fibers  and  matrix.  The  most  common  mode  for  practical  fiber  composites  is  the 
shear  mode.  Microbuckling  doesn’t  necessarily  damage  the  fibers  unless  it  is  followed 
by  kinking  of  the  fibers.  A  schematic  of  fiber  kinking  and  subsequent  fracture  is 
shown  in  Figure  3.  Kink  band  formation  has  been  observed  in  Kevlar,  glass,  and 
graphite  composites,  as  discussed  in  References  [12-13]. 

Microbuckling  with  Unbonded  Fibers.  A  failure  mode  for  composites  with  a 
weak  fiber-matrix  interface  is  fiber-matrix  debonding  followed  by  microbuckling  in 
the  shear  mode.  Microbuckling  theory,  and  tests  conducted  by  Greszczuk  on  ideal 
composites  (metal  rods  in  a  matrix)  indicate  that  the  strength  of  composites  with 
unbonded  fibers  is  substantially  less  than  that  with  bonded  fibers.  Hahn  et  al.  in 
Reference  [12]  also  observed  that  microbuckling  is  initiated  at  regions  of  local  inho- 
mogenities,  such  as  voids,  stress  concentrations,  or  weak  matrix  areas.  .Microbuckling 
of  unbonded  fibers  can  also  be  followed  by  fiber  kinking. 
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EXTENSIONAL  MODE  SHEAR  MODE 

Figure  2.  Microbuckling  failure  modes. 


Figure  3.  Kink  band  formation  and  fracture  following  microbuckling 
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Figure  4.  Interaction  failure  mechanism  in  fiber  composites. 

Interaction  Failure.  When  a  unidirectional  composite  is  loaded  in  longitudinal 
compression,  transverse  tensile  and  compressive  stresses  develop  in  the  matrix  because 
of  the  differences  in  Poisson’s  ratios  between  the  fibers  and  the  matrix,  as  shown 
in  Figure  4.  Consider  a  compressive  load  applied  to  unidirectional  composite  in 
the  fiber  direction.  Stresses  develop  in  the  matrix  in  the  transverse-fiber  directions. 
The  distribution  of  stress  in  the  matrix  adjacent  to  a  single  fiber,  called  a  repeating 
element,  is  shown  in  Figure  4  for  the  case  when  Poisson’s  ratio  of  the  fiber,  i/f,  is 
less  than  Poisson’s  ratio  of  the  matrix,  If  the  matrix  transverse  tensile  strength 
is  low,  matrix  cracking,  fiber-matrix  debonding,  or  delaminations  may  occur.  .As  an 
example,  test  data  are  compared  in  Figure  5  with  a  theory  developed  by  Greszczuk 
and  Pagano.  The  theory  relates  the  ratio  of  the  ultimate  compression  strength  with 
interaction,  Flc,  and  compression  strength  without  interaction,  F{kf  (Ff  is  the  fiber 
compressive  strength  and  k(  is  the  fiber  volume  fraction),  to  the  transverse  tensile 
strength,  Fjt. 

Direct  Fiber  Failure.  Fiber  compression  failure  occurs  before  microbuckling  if  the 
fibers  are  weak  in  compression. 


.A  comparison  of  test  data  with  calculated  stresses  for  each  of  the  compressive 
failure  modes  just  discus.sed  is  shown  in  Table  1.  This  table  is  reproduced  from 
Reference  [6].  Note  that  the  measured  failure  strength  of  Kevlar-epoxy  is  close  to 
that  calculated  from  microbuckling  theory  for  unbonded  fibers.  Fiber  kinking  follow¬ 
ing  microbuckling  is  also  reported  for  Kevlar-epoxy  in  Reference  [13].  .A  flow  chart 


Figure  5.  Comparison  of  interaction  failure  theory  with  test  data. 


(reproduced  from  Reference  [12])  which  describes  the  microstructural  compression 
response  of  fiber  composites  is  shown  in  Figure  6.  The  failure  sequence  suggested  for 
unidirectional  graphite  and  glass  reinforced  composites  is  given  by  the  shaded  boxes. 
The  chart  suggests  that  microbuckling,  followed  by  fiber  kinking  or  bending  failure 
due  to  continued  microbuckling,  are  the  dominate  compressive  failure  modes  for  these 
composites. 


1.1.3  Matrix  Failure  due  to  Transverse  Fiber  Stress  and 
Shear  Stress. 

The  primary  factors  affecting  the  transverse  tensile,  compressive,  and  shear  strengths 
of  a  composite  are: 

•  properties  of  the  fiber  and  matrix  (such  as  transverse  Young’s  moduli  for  trans¬ 
verse  failure  and  shear  moduli  for  shear  failure), 

•  fiber  volume  fraction, 

•  void  content, 

•  fiber-matrix  bond  strength, 

•  stress  concentrations  from  bonded  fibers  and  from  voids. 
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Table  1 .  Comparison  of  compressive  stresses  from  four  theoretical  failure  modes  with 
data  from  composites  made  with  rigid  resins. 


Measured  Data  vs.  Theoretical  Modes 

Composites 

Unbonded 

Bonded 

Modmor  II 

Kevlar 

Steel  Rods 

Steel  Rods 

Epoxy 

Epoxy 

in  Epoxy 

in  Epoxv 

(MPa) 

(MPa) 

(.MPa) 

(MPa) 

Measured  Data: 

481 

887 

972 

276 

Theoretical  Failure  Modes: 

Microbuckling  (Bonded  Fibers) 

— 

6255 

4552 

2069 

Microbuckling  (Unbonded  Fibers) 

533 

~ 

372 

372 

Interaction  Failure 

— 

881 

1200 

1207 

Direct  Fiber  Failure 

— 

946 

1338 

1517 

Mlcrottructural  Compratdon  Ratponaa 


Figure  6.  Micromechanical  compression  failure  modes  of  fiber  composites.  The  shaded 
boxes  give  the  failure  sequence  suggested  for  unidirectional  graphite  and 
gla.ss  reinforced  composites. 
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TRANSVERSE  TENSILE  STRENGTH  OF 


Figure  7.  Comparison  of  test  data  with  theory  for  transverse  tensile  failure  of 
graphite-epoxy. 


Greszczuk  in  Reference  [6]  shows  that  the  interaction  of  triaxial  stress  concentra¬ 
tions  from  fibers  and  voids,  and  the  number  of  unbonded  or  ineffective  fibers  aligned 
at  the  boundary  of  voids,  have  a  strong  effect  on  transverse  tensile  strength.  As 
an  example,  the  transverse  tensile  strengths  of  composites  made  with  graphite  fibers 
and  different  resins  are  compared  in  Figure  7  with  a  theory  developed  by  Greszczuk. 
The  abscissa  parameter,  /?,  depends  on  the  volume  fractions  of  fibers  (kf).  voids  (kv). 
and  ineffective  fibers  (kjf).  Good  agreement  between  test  data  and  theory  is  also  re¬ 
ported  by  Greszczuk  for  S-glass  fibers  embedded  in  different  resins  (not  shown).  For 
S-glass  composites,  Greszczuk  shows  that  the  transverse  tensile  strength  approaches 
(and  exceeds  slightly)  the  tensile  strength  of  the  resin  as  /?  — +  0,  i.t.  cis  the  void  and 
ineffective-fiber  fractions  approach  zero. 

Tensile  tests  in  the  transverse  fiber  direction  of  unidirectional  specimens  indicate 
that  cracks  form  in  the  matrix  parallel  to  the  fiber  direction,  as  discussed  in  Refer¬ 
ence  [8]  and  shown  schematically  in  Figure  8.  Both  transverse  fiber  and  shear  stresses 
are  assumed  to  contribute  to  the  tensile  matrix  failure  mode.  For  transverse  fiber 
compressive  stress,  matrix  compression  failure  occurs.  Compression  failure  occurs  by 
shearing  in  the  matrix  in  a  plane  parallel  to  the  fibers,  but  not  necessarily  normal  to 
the  direction  of  the  applied  compressive  stress.  A  schematic  of  matrix  compression 
failure  is  shown  in  Figure  9.  Both  transverse  fiber  compressive  stress  and  shear  stress 
are  assumed  to  contribute  to  the  matrix  compressive  failure  mode. 
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Figure  8.  Matrix  tensile  failure  mode  in  unidirectional  composite. 


Figure  9.  Matrix  compression  failure  mode  in  unidirectional  composite 


1.1.4  Comments  on  Strength  Coupling. 

As  discussed  in  the  preceeding  sections,  both  tensile  and  compressive  fiber  failure 
strengths  are  coupled  with  the  matrix  and  fiber-matrix  debond  strengths.  Matrix 
failure  reduces  the  ability  of  the  matrix  to  transfer  load  through  the  lamina,  which 
reduces  the  fiber  tensile  failure  strength  of  the  composite.  Fiber  composites  can  also 
fail  in  longitudinal  compression  due  to  low  matrix  strength.  Strength  coupling  is  also 
path  dependent.  For  example,  Wu  in  Reference  [6]  suggests  that  broken  fibers  can 
initiate  transverse  matrix  cracks,  thus  reducing  the  transverse  tensile  strength  of  the 
composite.  Biaxial  stress  tests  should  be  performed  on  unidirectional  composites  to 
aid  in  understanding  strength  coupling  and  to  validate  interactive  failure  criteria. 

1.2  REVIEW  OF  FAILURE  CRITERIA. 

Failure  criteria  are  used  to  relate  critical  combinations  of  stresses  or  strains  to  fail¬ 
ure  in  a  material.  Numerous  failure  criteria  are  available  to  predict  the  strength  of 
laminated  materials.  References  [14-20].  Simple  criteria  are  applied  directly  to  the 
total  laminate,  but  more  complete  criteria  are  applied  to  the  individual  layers  of  the 
laminate. 

In  direct  laminate  failure  criteria,  the  laminate  is  considered  to  be  homogeneous 
and  anisotropic.  Laminatidn  theory  is  not  required  to  analyze  the  stresses  and  strains 
in  each  layer.  However,  the  criteria  are  difficult  to  use  because  laminate  strengths 
must  be  determined  for  each  layup  that  is  analyzed. 

Lamina-by- lamina  failure  criteria  are  more  common  than  laminate  failure  criteria 
because  they  are  more  general,  and  hence  more  universally  applicable.  Only  one  set  of 
material  property  tests  need  be  performed  on  unidirectional  specimens  to  determine 
the  strength  characteristics  of  the  fiber  composite  lamina,  and  then  any  layup  can  be 
analyzed.  Lamination  theory  is  used  to  obtain  the  stresses  and  strains  in  the  principal 
material  directions  of  each  lamina.  The  failure  criteria  are  applied  to  each  lamina  to 
predict  lamina-by-lamina  failure. 

Lamina  failure  criteria  are  macromechanical  because  the  average  stress  and  strain 
in  the  lamina  are  used  to  predict  failure.  The  lamina  are  considered  to  be  homo¬ 
geneous  and  orthotropic.  Micromechanical  failure  criteria  are  less  common  than 
macromechanical  failure  criteria  because  of  the  complexity  of  dealing  with  the  inter¬ 
action  of  the  fiber  and  matrix  components.  However,  guidelines  based  on  microme¬ 
chanics,  such  35  understanding  how  the  stresses  in  the  fiber,  matrix,  and  fiber-matrix 
interface  each  contribute  to  failure,  are  used  to  establish  failure  criteria  that  distin¬ 
guish  between  failure  modes. 

Sometimes  micromechanics  is  used  to  predict  the  strength  of  the  lamina  from  the 
strength  of  the  components,  but  more  often,  lamina  strengths  are  obtained  experi¬ 
mentally  from  uniaxial  tests  on  unidirectional  specimens.  A  few  criteria  also  require 
biaxial  test  data.  The  strengths  measured  for  use  in  a  plane-stress  failure  criteria 
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are  generally  the  tensile  and  compressive  strengths  in  the  fiber  direction,  the  tensile 
and  compressive  strengths  in  the  transverse  fiber  direction,  and  the  in-plane  shear 
strength.  Failure  modes  of  unidirectional  fiber-composites  are  discussed  in  Section 
1.2. 

The  failure  surface  of  an  orthotropic  lamina  in  a  state  of  plane  stress  is  three- 
dimensional.  It  is  defined  in  stress  space  by  the  stresses  in  the  principal  material 
directions  of  the  lamina,  an,  <T22,  and  aij  (ai,  a2,  and  ae  in  contracted  notation).  The 
intercepts  of  the  failure  surface  with  the  coordinate  axes  are  the  material  strengths 
measured  from  uniaxial  tests  on  unidirectional  specimens.  For  simplicity,  the  failure 
envelope  can  be  plotted  as  a  curve  in  Cn  ,  a22  stress  space,  for  constant  values  of 
(Ti2-  Most  criteria  are  in  agreement  for  the  failure  stresses  in  the  principal  material 
directions,  i.e.  the  intercepts  of  the  failure  surfaces  with  the  coordinate  axes  are  in 
agreement.  The  criteria  differ  on  what  constitutes  failure  for  biaxial  stress  states. 

A  survey  of  the  most  commonly  used  failure  criteria  was  conducted  by  the  Amer¬ 
ican  Institute  of  Aeronautics  and  Astronautics  (AIAA)  to  determine  how  much  vari¬ 
ability  exists  in  the  design  and  analysis  of  composite  structures,  Reference  [14].  The 
survey  results  indicate  that  the  maximum  stress,  maximum  strain,  Tsai-Hill,  and 
Tsai-Wu  failure  criteria  are  most  often  used  by  surveyed  AIAA  members.  These 
failure  criteria  are  discussed  in  the  following  paragraphs. 

Maximum  Stress.  Failure  occurs  when  any  component  of  stress  in  the  principal 
material  directions  of  the  lamina  exceeds  its  corresponding  strength: 
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The  five  ultimate  strengths  are  determined  from  uniaxial  and  pure-shear  tests  on 
unidirectional  specimens: 

Xt  Tensile  strength  in  longitudinal  fiber  direction 
Xc  Compressive  strength  in  longitudinal  fiber  direction 
Yt  Tensile  strength  in  transverse  fiber  direction 
Yc  Compressive  strength  in  transverse  fiber  direction 
5  In-plane  shear  strength 

Five  independent  modes  of  failure  are  predicted:  tensile  and  compressive  longitudinal 
failure,  tensile  and  compressive  transverse  failure,  and  shear  failure. 

Maximum  Strain.  Failure  occurs  when  any  strain  component  in  the  principal  ma¬ 
terial  directions  of  the  lamina  exceeds  its  corresponding  ultimate  strain.  Five  ultimate 
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strain  values  are  determined  from  uniaxial  and  pure  shear  tests  on  unidirectional  spec¬ 
imens;  the  ultimate  tensile  and  compressive  strains  in  the  fiber  and  transverse  fiber 
directions,  and  the  ultimate  shear  strain.  Five  failure  modes  are  also  predicted:  ten¬ 
sile  and  compressive  failure  in  the  longitudinal  fiber  and  transverse  fiber  directions, 
and  shear  failure.  These  are  the  same  failure  modes  predicted  with  the  maximum 
stress  criteria.  The  maximum  strain  failure  criteria  allow  for  some  interaction  among 
the  in-planj  stresses  through  the  Poisson  effect. 

Tsai-Hill.  This  failure  criterion  is  a  generalization  of  the  maximum  distortional 
energy  criterion  for  isotropic  materials.  Hill,  in  Reference  [15],  extended  the  criterion 
to  orthotropic  materials.  Tsai  and  Azzi,  in  Reference  [16],  modified  the  Hill  criterion 
to  account  for  the  plane  stress  conditions  and  transverse  isotropy  of  fiber  composite 
lamina: 

<5‘u  <^11  <^22  <^22  ,  ^12  ^ 

J^2  Y2  y^2  _5'2  —  ^  ^  ' 

where  X  and  Y  are  the  in-piane  strengths  in  the  fiber  and  transverse  fiber  directions, 
respectively,  and  S  is  the  in-plane  shear  strength.  Tsai  showed  that  the  criterion  is 
applicable  to  composites  with  different  properties  in  tension  and  compression.  Ten¬ 
sile  strengths  are  used  when  the  corresponding  stresses  are  tensile,  and  compressive 
strengths  are  used  when  the  corresponding  stresses  are  compressive.  Tsai  also  devel¬ 
oped  two  additional  equations  for  mutually  orthogonal  planes,  similar  to  Equation  ( 1 ). 
for  failure  analysis  of  three-dimensional  materials.  Reference  [17]. 

The  Tsai-Hill  criterion  is  called  an  interactive  failure  criteria  because  failure  de¬ 
pends  on  more  than  one  value  of  stress.  The  failure  surface  for  the  Tsai-Hill  criterion 
is  a  smooth  surface  (ellipsoidal  or  spherical,  depending  on  the  strength  values)  in 
f^22t  <^12  stress  space.  In  contrast,  the  failure  surfaces  of  the  maximum  stress  and 
maximum  strain  criteria  are  not  smooth  surfaces  because  they  are  formed  by  inter¬ 
secting  planes.  Unlike  the  maximum  stress  and  strain  criteria,  the  mode  of  failure 
is  not  predicted  by  the  Tsai-Hill  criterion.  These  two  types  of  criteria  are  analo¬ 
gous  to  Tresca-type  (intersecting  planes)  and  Mi^es-type  (smooth  surface)  criteria  for 
homogeneous  isotropic  materials. 

Tsai-Wu.  Tsai  and  Wu  developed  a  tensor-component  polynomial  theory  as  a  fail¬ 
ure  criterion  for  anisotropic  materials.  Reference  [18].  The  criterion  contains  linear  as 
well  as  quadratic  stress  terms.  Failure  occurs  when  the  following  equation  is  satisfied: 

FiCTii  -h  F2<722  -h  FiiCTjj  -|-  Fl|2Cr22  +  F\2C^u'^22  +  ^66(712  >1  (-) 

Six  coefficients  must  be  defined  for  a  composite  lamina  in  a  state  of  plane  stress. 
Five  coefficients  {Fi,  F2,  Fn,  F22,  Fee)  are  determined  from  uniaxial  and  shear  tests 
on  unidirectional  specimens.  Each  of  the  coefficients.  Fi,  F2,  Fu,  and  F22  include 
contributions  from  both  tensile  and  compressive  strengths  (see,  for  example.  Equa¬ 
tion  26),  rather  than  requiring  separate  values  for  tension  and  compression  as  for  the 
.V  and  Y  strengths  in  the  Tsai-Hill  criterion.  The  si.xth  coefficient,  F12,  must  be  de¬ 
termined  from  biaxial  tests,  a  variety  of  which  are  available.  Unfortunately,  different 
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(a)  Tsai-Wu  Polynomial  RMS=0.137 
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!b) Maximum  Strain  RMS=0.200 


02 
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Figure  10.  Comparison  of  three  commonly  used  failure  criteria  with  graphite-epoxy 
test  data. 

tests  produce  different  values  of  Fi2-  Like  the  Tsai-Hill  failure  criterion,  the  failure 
envelope  is  a  smooth  surface  in  stress  space  and  only  the  onset  of  failure  is  predicted, 
not  the  mode  of  failure. 


The  .4IAA  survey  indicated  that  the  maximum  strain  criterion  is  the  most  fre¬ 
quently  used  failure  criterion,  followed  by  the  maximum  stress,  Tsai-Hill.  and  Tsai-Wu 
criteria.  The  ability  of  each  criterion  to  predict  experimental  data  varies  from  one 
composite  material  to  another.  Comparison  of  three  of  the  failure  criteria  with  test 
data  for  graphite-epoxy  is  shown  in  Figure  10.  This  figure  is  reproduced  from  Ref¬ 
erence  [19].  The  root  mean  square  (RMS)  values  indicate  how  well  the  criteria  fit 
the  test  data.  An  advantage  of  the  maximum  strain  and  maximum  stress  criteria 
is  that  the  mode  of  failure  is  predicted.  With  the  mode  of  failure  identified,  pro¬ 
gressive  damage  in  the  laminate  can  be  analyzed  by  including  degradation  of  lamina 
properties. 

The  Chang  criteria,  which  are  currently  implemented  in  DYN.A.3D.  are  not  com¬ 
pared  with  test  data  in  Figure  10,  because  material  strengths  and  failure  modes  were 
not  reported  for  the  test  data.  The  Chang  criteria  are  discussed  in  detail  in  Section 
1.3. 
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1.3  DEVELOPMENT  OF  CHANG’S  FAILURE 
CRITERIA. 

As  discussed  in  Section  1.1,  the  Tsai- Hill  and  Tsai-Wu  interactive  failure  criteria  pre¬ 
dict  when  a  given  set  of  stresses  and  strains  will  produce  failure,  but  they  do  not 
give  the  mode  of  failure.  Hashin  developed  an  interactive  failure  criterion,  Refer¬ 
ence  [2],  in  which  four  distinct  failure  modes  are  modeled;  tensile  and  compressive 
fiber  failure  and  tensile  and  compressive  matrix  failure.  Chang  modified  Hashin's 
criterion  to  include  nonlinear  shear  stress-strain  behavior.  Chang  also  defined  a  post¬ 
failure  degradation  rule  so  that  the  behavior  of  the  laminate  can  be  analyzed  after 
each  sucessive  lamina  fails.  The  Chang  failure  criterion.  References  [3-5]  is  used  in 
the  DYNA3D  composite  damage  model.  In  the  following  discussion,  Hashin’s  fail¬ 
ure  criterion,  and  then  Chang’s  modification  to  include  nonlinear  shear  behavior,  are 
presented.  Finally,  the  DYNA3D  post-failure  degradation  rule  is  summarized. 

1.3.1  Hashin  Failure  Criteria. 

Hashin’s  three-dimensional  failure  criteria  were  established  for  unidirectional  fiber 
composites  in  terms  of  a  quadratic  stress  polynomial.  Since  most  fiber  composites 
are  transversely  isotropic,  Hashin  first  defined  a  general  failure  criterion  in  terms 
of  the  stress  invariants  of  a  transversely  isotropic  material.  The  three-dimensional 
criterion  has  the  general  form: 

A\Ii  -f  Bxl^  -I-  A2I2  +  B2I2  +  C'12/1/2  +  A3/3  -f-  A4I4  >  1  (3) 

where  /i,  /2,  /s,  and  are  the  stress  invariants; 

A  =  cr,i 

A  =  <^22  +  <^33 

A  =  —  Cr22Cr33  (4) 

A  =  <^12  d"  ^13 

Hashin  simplified  the  criterion  for  fiber  composites  in  a  state  of  plane  stress. 

Hashin  argued  that  failure  is  produced  by  the  normal  and  shear  stresses  acting  on  the 
failure  plane.  For  failure,  the  failure  plane  is  the  2-3  plane,  acted  on  by  stresses 
(7ii  and  <ti2.  The  stress  in  the  transverse  fiber  direction  [(722)  does  not  contribute  to 
fiber  failure.  In  matrix  failure,  a  plane  crack  occurs  parallel  to  the  fiber  direction. 
Failure  occurs  in  any  plane  with  axes  parallel  and  transverse  to  the  fibers.  The  matrix 
failure  plane  is  acted  on  by  stresses  <722  and  ctjj.  He  argued  that  the  stress  in  the 
longitudinal  fiber  direction  (<Tn)  does  not  contribute  to  matrix  failure  because  this 
stress  is  carried  almost  entirely  by  the  fibers. 

Hashin’s  plane  stress  failure  criteria,  defined  in  terms  of  the  five  principal  strengths 
of  a  fiber  composite,  are: 
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Fiber  Tensile  Failure:  Hashin  assumed  that  tensile  and  shear  stresses  are  mutu¬ 
ally  weakening,  hence  both  contribute  to  fiber  breakage.  This  criterion  is  also  known 
as  the  Yamada-Sun  failure  criterion,  Reference  [21].  The  Yamada-Sun  failure  crite¬ 
rion  is  a  simplified  form  of  the  Tsai-Hill  criterion.  Equation  1,  expressed  with  the 
transverse  fiber  stress  set  equal  to  zero  (<r22  =  0). 

©  ■^(t) 

Fiber  Compressive  Failure:  Hashin  represents  fiber  compressive  failure  in  simple 
maximum  stress  form.  Hashin  argued  that  there  is  no  physically  reasonable  method 
for  including  the  effect  of  shear  stress. 


(Tn  >  Xc  (fTii  <  0) 


(6) 


Matrix  Tensile  Failure:  Hashin’s  matrix  tensile  failure  criterion  is  a  simplified 
form  of  the  Tsai-Hill  criterion.  Equation  1,  expressed  with  the  longitudinal  fiber 
stress  set  equal  to  zero  (<Tu  =0). 


Matrix  Compressive  Failure:  Matrix  failure  in  the  transverse  fiber  direction 
depends  on  the  compressive  strength  in  the  transverse  fiber  direction,  Yc,  the  in¬ 
plane  shear  strength,  5,  and  the  transverse  shear  strength,  St-  Hashin’s  matrix 
compression  failure  criterion  contains  a  first-order  term  in  <722  whose  sign  depends  on 
Yc  and  Sj-  The  effect  of  this  linear  term  on  the  lamina  failure  envelope  is  discussed 
in  Section  2.2. 


((722  <  0) 


(S) 


1.3.2  Chang  Failure  Criteria  used  in  DYNA3D  Composite 
Damage  Model. 


Chang  et  al. ,  References  [3-5],  used  a  modified  form  of  Hashin 's  fiber  breakage,  ma¬ 
trix  tensile,  and  matrix  compressive  failure  criteria  to  analyze  the  failure  strength  of 
laminates  exhibiting  nonlinear  shear  behavior.  Chang  modified  each  of  the  failure 
criteria  to  include  nonlinear  shear  effects,  following  Sandhu's  strain  energy  method. 
Reference  [22].  The  modified  form  of  the  fiber  breakage  failure  criterion  (Equation  5) 
is; 


+ 


n\2 

(7i2cf7i2 

T"  <^l2dlu 

Jo 


>  1 


(9) 
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where  712  is  the  ultimate  shear  strain  of  the  lamina.  Chang  used  the  nonlinear  shear 
stress-strain  relation  advanced  by  Hahn  and  Tsai,  Reference  [23],  in  Equation  9: 


1 

7i2  =  7r-<^i2  +  OICT 
(J12 


3 

12 


(10) 


where  a  is  the  nonlinear  shear  stress  parameter.  This  parameter  is  included  in  the 
sheaj  stress-strain  formula  to  describe  the  nonlinear  behavior  typical  of  many  fiber- 
composite  materials.  The  parameter  a  is  determined  from  materials  testing.  We  refer 
to  the  term  containing  the  shear  stress-strain  integrals  as  the  fiber-matrix  shearing 
term.  By  performing  the  indicated  integration,  the  fiber-matrix  shearing  term,  f,  is 
defined  cis  follows; 


\  S  )  \  l  +  laGnS^  J 


(11) 


For  a  =  0,  this  term  reduces  to  the  ratio  of  the  shear  stress  to  shear  strength,  [ctu/ Sy, 
used  in  the  Hashin  criteria.  Equations  5,  7,  and  8.  For  o  >  0,  the  term  in  brackets  is 
less  than  one  so  that  f  <  {(Tx^lSy.  Hashin’s  tensile  and  compressive  matrix  failure 
models  are  also  modified  in  a  similar  manner  to  account  for  nonlinear  shearing. 


When  a  failure  criterion  is  met,  some  or  all  of  the  lamina’s  elastic  constants  are 
set  to  zero,  following  a  method  Chang  pursued  in  early  unpublished  works.  These 
degradation  rules  were  made  available  to  us  through  personal  communication  with 
Dr.  Chang.  More  recently,  Chang  improved  on  his  post-failure  degradation  rules  and 
published  results  in  References  [3-5].  We  will  therefore  refer  to  the  degradation  rules 
implemented  in  DYNA3D  as  the  DYNA3D  post-failure  model. 

The  DYNA3D  post-failure  model  is  an  ad  hoc  treatment  of  post-failure  degrada¬ 
tion  and  is  characterized  as  a  ‘strain  softening’  constitutive  response,  e.g.  after  failure 
the  strain  continues  to  increase  under  constant  or  decreasing  load.  One  fundamen¬ 
tal  problem  with  the  implementation  of  strain  softening  models  into  finite  element 
codes  is  mesh  size  dependency,  I'.e.  the  same  physical  problem  will  produce  differ¬ 
ent  results  for  different  mesh  configurations.  Strain  softening  effects  in  the  current 
composite  damage  model  version  of  DYNA3D  have  not  been  investigated.  However, 
Simo  in  Reference  [24]  developed  a  mathematically  rigorous  approach  to  strain  soft¬ 
ening,  which  he  calls  the  principal  of  Maximum  Damage  Dissipation.  His  work  was 
supported,  in  part,  by  a  subcontract  from  APTEK  and  may  be  utilized  in  future 
efforts. 


The  complete  failure  criteria  advanced  by  Chang  and  DYNA3D  post-failure  degra¬ 
dation  rules  are  defined  in  the  following  paragraphs.  In  his  criteria,  Chang  suggests 
using  the  in  situ  in-plane  shear  strength  measured  from  tests  on  symmetric  cross- 
ply  laminates,  rather  than  the  shear  strength  measured  from  tests  on  unidirectional 
specimens.  Chang  and  Springer  in  Reference  [4]  state  that  the  shear  strength  of  a 
three-layered  glass-epoxy  composite  is  two  to  three  times  greater  than  that  obtained 
by  testing  single  plies. 
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(12) 


Fiber  Breakage  Failure  Criterion; 

When  the  fiber  breakage  failure  criterion  is  satisfied,  all  of  the  elastic  constants  of  the 
failed  lamina  are  set  to  zero,  i.e.  Ell  —  E22  =  Gi2  =  t'i2  =  t/21  =  0.  The  lamina  is  no 
longer  able  to  carry  load  in  the  fiber,  transverse  fiber,  or  shear  directions.  The  load 
is  redistributed  to  the  remaining  lamina. 

Matrix  Cracking  Failure  Criterion: 

(^)'  +  f>l  ((722  >0)  (13) 

When  the  matrix  cracking  failure  criterion  is  satisfied,  all  of  the  elastic  constants 
except  for  the  fiber  modulus,  Fu,  are  set  to  zero,  i.e.  E22  =  G12  =  ^^12  =  ^21  =  0. 
The  lamina  is  unable  to  carry  load  in  the  transverse  fiber  or  shear  directions.  Once 
matrix  cracking  failure  occurs,  both  the  stress  in  the  transverse  fiber  direction  and 
shear  stress  are  reduced  to  zero. 


Matrix  Compression  Failure  Criterion: 


V25>' 


+ 


(i)'- 


(722 


+  r  >  1  ((722  <  0) 


(14) 


When  the  matrix  compression  failure  criterion  is  satisfied,  the  transverse  elastic  mod¬ 
ulus  and  both  Poisson’s  ratios  are  set  to  zero,  i.e.  E22  =  1^12  =  ^21  =  0-  The  lamina 
is  unable  to  carry  load  in  the  transverse  fiber  direction;  the  stress  in  the  transverse 
fiber  direction  is  reduced  to  zero.  Note  that  in  Equation  14  Chang  has  replaced  the 
transverse  shear  strength,  5r,  in  Kashin’s  criterion.  Equation  8,  with  the  in-plane 
shear  strength,  S. 


The  above  failure  criteria  and  post-failure  degradation  rules,  which  we  refer  to 
as  the  Chang  failure  criteria,  are  used  in  the  DYNA3D  composite  damage  model. 
The  Chang  criteria,  along  with  the  other  fiber-composite  failure  criteria^  are  given 
in  Table  2.  The  Chang  criteria  do  not  model  all  types  of  damage  observed  in  fiber 
composite  materials,  e.g.  fiber  compression  failure  and  delamination.  Static  fiber- 
direction  compression  tests  on  Kevlar-epoxy  indicate  that  the  stress-strain  response 
curves  can  be  modeled  cis  elastic-perfectly  plastic,  as  suggested  in  References  [25] 
and  [26].  The  probable  failure  mode  is  microbuckling  of  debonded  fibers.  .A  fiber 
compressive  failure  model  will  be  implemented  in  APTEK’s  version  of  DYN.4.3D  once 
an  appropriate  degradation  rule  is  established.  Dr.  Hallquist  has  also  implemented 
into  DYNA3D  a  failure  criterion  ba^ed  on  a  free-edge  delamination  model  for  the 
composite  brick  elements,  but  the  criterion  has  not  been  verified  by  .APTEK. 

'The  modified  Tsai-Wu  criterion  for  matrix  tensile  and  compressive  failure  is  discussed  in  Section 
2..3.3. 


Table  2.  Failure  criteria  for  fiber-composite  materials. 


Maximum  Stress: 
Fiber  Tension 

Fiber  Compression 

Matrix  Tension 

Matrix  Compression 

Matrix  Shear 
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Modified  Tsai-Wu: 
Matrix  Failure 


Hashin: 

Fiber  Tension 

Fiber  Compression 
Matrix  Tension 
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Chang: 

Fiber  Breakage 

Matrix  Cracking 
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SECTION  2 

VERIFICATION  OF  DYNA3D  COMPOSITE 
SHELL  DAMAGE  MODEL. 


Calculations  were  performed  to  verify  the  implementation  of  the  failure  criteria  and 
degradation  rules  used  in  the  DYNA3D  composite  damage  model.  Each  failure  mode 
and  degradation  rule  is  demonstrated  through  calculations  on  single-element  shell 
laminates.  Computational  results,  using  the  damage  model,  are  also  compared  with 
experimental  results  found  in  the  literature  to  verify  the  theoretical  basis  of  the  failure 
criteria  and  post-failure  degradation  rules.  Theoretical  aspects  of  the  failure  criteria 
were  also  evaluated  by  Dr.  Juan  Simo  and  Dr.  Junho  Jang  of  Stanford  University, 
under  subcontract  to  APTEK.  Their  efforts  are  included  in  two  separate  reports. 
References  [27]  and  [24]. 

2.1  CALCULATIONS  TO  DEMONSTRATE  FAILURE 
MODES. 

The  Chang  composite  damage  model  in  DYNA3D  treats  three  modes  of  damage: 
fiber  breakage  (tensile  failure),  matrix  cracking  (tensile  failure),  and  matrix  compres¬ 
sion  failure.  Four  material  strength  parameters  are  used  in  the  three  failure  criteria; 
longitudinal  tensile  strength  A'x,  transverse  tensile  strength  Yr,  transverse  compres¬ 
sive  strength  Yc,  and  shear  strength  S.  A  fifth  parameter,  the  nonlinear  shear  stress 
parameter,  a,  is  also  included  in  the  damage  model  to  describe  the  nonlinear  in-plane 
shear  stress-strain  behavior.  For  each  damage  mode,  failure  in  a  lamina  occurs  when  a 
combination  of  the  stresses  exceeds  the  appropriate  strength  relation  for  the  material. 
A  mathematical  description  of  the  damage  model  was  included  in  Section  1.3. 

Calculations  were  performed  to  verify  the  implementation  of  the  fiber  and  ma¬ 
trix  failure  modes  in  DYNA3D.  Velocities  were  applied  to  each  node  of  an  element 
composed  of  a  single  orthotropic  lamina  to  achieve  a  state  of  combined  extensional 
and  shear  strain.  The  strains  in  the  element  increased  steadily  until  failure  occured. 
The  material  strength  parameters  were  selected  so  that  lamina  failure  in  each  of  the 
damage  modes  would  occur  at  a  specified  strain  level,  corresponding  to  a  known  stress 
level.  The  purpose  of  these  calculations  Wcis  to  verify  the  implementation  for  each 
of  the  three  failure  criteria  by  comparing  the  stress  level  at  failure  with  the  com¬ 
puted  failure  stress;  we  are  satisfied  that  the  failure  criteria  in  our  DYN.\3D  version 
are  implemented  correctly.  Elastic  constants  and  strengths  used  in  the  verification 
calculations  are  listed  in  Table  3. 
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Table  3.  Material  properties  used  in  verification  calculations. 


£^11 

= 

51,710 

MPa 

(7500  ksi) 

E22 

= 

6,897 

MPa 

(1000  ksi) 

G\2 

=: 

3,448 

MPa 

(  500  ksi) 

V\2 

= 

0.3 

Xt 

= 

1379 

Mpa 

(200  ksi) 

Yt 

= 

345 

Mpa 

(  50  ksi) 

Vc 

= 

69 

Mpa 

(  10  ksi) 

s 

= 

83 

Mpa 

(  12  ksi) 

2.1.1  Fiber  Failure  Modes. 

Fiber  Breakage.  For  fiber  breakage  to  occur,  the  stress  in  the  fiber  direction  of 
the  lamina  must  be  tensile.  When  the  fiber  breakage  failure  criterion  is  satisfied, 
all  of  the  in-plane  elastic  constants  for  that  particular  lamina  are  degraded,  i.e.  E\\, 
£•221  G^i2i  and  t'ji-  The  elastic  constants  are  degraded  over  an  arbitrary  number 
of  time  steps  (100  time  steps  were  used  in  this  demonstration)  until  they  are  equal 
to  zero.  Gradual  degradation  of  the  elastic  constants  was  implemented  to  minimize 
dynamic  instabilities  associated  with  sudden  load  redistribution. 

Results  from  a  calculation  demonstrating  fiber  breakage  are  shown  in  Figure  11. 
The  calculation  was  performed  on  a  single-element  shell  laminate  with  unidirectional 
fibers  (oriented  in  the  x-direction).  Monotonically  increasing  shear  and  longitudinal 
(in  the  direction  of  the  fibers)  strains  were  applied  to  the  laminate  simultaneously. 
The  ratio  of  the  applied  shear  strain  to  the  applied  longitudinal  strain  wa^  one.  Fiber 
breakage  occurs  when  a  combination  of  the  stresses  exceeds  the  fiber  breakage  strength 
relation,  Equation  12.  The  failure  stresses  are  <7n  =  825  MPa  and  =  "1  MPa. 
After  failure,  both  the  stress  in  the  fiber  direction  and  the  shear  stress  are  reduced 
to  zero  over  100  time  steps  while  the  corresponding  strains  increaise.  Note  that  the 
nonlinearity  observed  in  the  shear  stress-strain  behavior  is  due  to  the  non-zero  shear 
stress  parameter  used  in  the  calculation  (a  =  3.05  x  10"*  MPa"*). 

Fiber  Compressive  Failure.  .At  present,  the  damage  model  does  not  check  the 
compressive  stress  in  the  fiber  direction  to  determine  if  fiber  compressive  failure  would 
occur.  To  include  such  failure,  one  would  have  to  determine  whether  or  not  the  failure 
criterion  and  post-failure  degradation  rule  depend  on  the  mode  of  compressive  failure. 
For  instance,  interactive  failure  may  affect  property  degradation  in  the  transverse 
fiber  and  shear  directions,  while  fiber  kinking  and  subsequent  fracture  may  also  affect 
Young’s  modulus  in  the  longitudinal-fiber  direction. 
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Extsnslonal  Strain  (percent) 

(a)  Stress-strain  behavior  in  longitudinal  fiber  direction. 


2.1.2  Matrix  Failure  Modes. 

Matrix  Cracking  Failure.  For  matrix  cracking  failure  to  occur,  the  stress  trans¬ 
verse  to  the  fiber  direction  must  be  tensile.  If  the  failure  criterion  is  satisfied,  the 
transverse  Young’s  modulus,  shear  modulus,  and  both  the  Poisson’s  ratios  are  set 
to  zero  for  that  particular  lamina,  i.e.  E22  =  Gu  =  1^12  =  ^^21  =  0.  At  present, 
these  moduh  are  degraded  over  a  single  time  step.  A  literature  search  was  performed 
to  determine  what  degradation  rules  are  appropriate  for  the  matrix  failure  modes. 
i.e.  which  elastic  constants  should  be  degraded.  Results  of  the  literature  search  and 
modifications  to  the  DYNA3D  implementation  are  discussed  in  Section  2.4. 

Results  from  a  calculation  demonstrating  matrix  cracking  failure  are  shown  in 
Figure  12.  The  calculation  was  performed  on  a  single-element  shell  laminate  with 
unidirection<d  fibers  (oriented  in  the  x-direction).  A  monotonically  increasing  strain 
was  applied  in  the  transverse  fiber  direction  (y  direction).  A  monotonically  increasing 
shear  strain  was  simultaneously  applied  to  the  launinate.  The  ratio  of  the  applied  shear 
strain  to  applied  longitudinal  strain  was  one.  A  nonlinear  shear  stress  parameter  of 
a  =  3.05  X  10~*  MPa"^  was  used  in  the  calculation.  When  a  combination  of  the 
transverse  stress  and  shear  stress  exceeds  the  strength  relation.  Equation  13.  matrix 
cracking  failure  occurs.  The  failure  stresses  are  <^22  =  130  Mpa  and  a-i2  =  78  MPa. 
Both  the  stress  in  the  transverse  fiber  direction  and  the  shear  stress  drop  to  zero  over 
one  time  step,  while  the  corresponding  strains  become  large. 

Matrix  Compression  Failure.  For  matrix  compression  failure  to  occur,  the  stress 
transverse  to  the  fiber  direction  must  be  compressive.  If  the  failure  criterion  is  satis¬ 
fied,  the  transverse  Young’s  modulus  and  both  the  Poisson’s  ratios  are  set  to  zero  for 
that  particular  lamina,  i.e.  E22  =  1^12  =  1^21  =  0- 

Results  from  a  calculation  demonstrating  matrix  compression  failure  are  shown 
in  Figure  13.  The  calculation  was  performed  on  a  single-element-shell  laminate  with 
unidirectional  fibers  (oriented  in  the  x-direction).  A  monotonically  increasing  com¬ 
pressive  strain  was  applied  in  the  transverse  fiber  direction  (y  direction).  A  monoton¬ 
ically  increasing  shear  strain  was  simultaneously  applied  to  the  laminate.  The  ratio  of 
the  applied  shear  strain  to  applied  longitudinal  strain  was  minus  one.  When  a  com¬ 
bination  of  the  transverse  compressive  stress  and  shear  stress  exceeds  the  strength 
relation.  Equation  14,  matrix  compression  failure  occurs.  The  failure  stresses  are 
<722  =  48  Mpa  and  0-12  =  48  MPa.  The  stress  in  the  transverse  fiber  direction  drops  to 
zero  over  one  time  step  while  the  corresponding  strain  increases.  However,  the  shear 
stress  continues  to  increase  with  increasing  shear  strain  because  the  shear  modulus 
is  not  degraded  in  the  Chang  post-failure  degradation  rule  for  matrix  compression 
failure. 
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(a)  Stress-strain  behavior  in  transverse  fiber  direction. 


Figure  12.  Demonstration  of  matrix  cracking  failure. 
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Extansional  Strain  (parcani) 

(a)  Stress-strain  behavior  in  transverse  fiber  direction. 


Shaar  Strain  (parcant) 

(b)  Shear  stress-strain  behavior. 


Figure  13.  Demonstration  of  matrix  compression  failure. 


2.2  CALCULATIONS  TO  DEMONSTRATE 
POST-FAILURE  BEHAVIOR. 

Two  calculations  were  performed  to  verify  the  post-failure  behavior  of  cross-ply  lam¬ 
inates  after  initial  ply  failure.  In  the  first  calculation  the  unloading  behavior  of  a 
cross-ply  laminate  was  examined  using  the  post-failure  degradation  rules  implemented 
in  APTEK’s  version  of  DYNA3D^;  we  are  satisfied  that  the  post-failure  degradation 
rules  are  implemented  correctly.  In  the  second  calculation  the  behavior  of  a  cross-ply 
laminate  which  exercises  a  different  post-failure  degradation  rule,  viz.  constant  stress 
after  matrix  failure,  was  examined.  Comparison  of  the  two  calculations  indicates 
that  the  ultimate  strength  of  the  laminate  is  extremely  sensitive  to  the  post-failure 
behavior  assumed  for  the  matrix. 

In  these  two  verification  calculations,  velocities  are  applied  to  each  node  of  an 
element  composed  of  alternating  cross-ply  lamina  to  achieve  a  state  of  combined 
extensional  and  shear  strain.  The  applied  extensional  strain  is  in  the  fiber  direction  of 
some  of  the  cross  plies  and  transverse  to  the  fiber  direction  of  the  remaining  cross  plies. 
The  material  strength  parameters  were  chosen  to  demonstrate  failure  and  degradation 
in  each  of  the  tensile  failure  modes.  Failure  and  material  property  degradation  depend 
on  the  fiber  orientation  of  each  lamina  and  the  sign  of  the  stress. 

2.2.1  Post-Failure  Behavior  of  a  Symmetric  Cross-Ply 
Laminate. 

Results  from  a  calculation  demonstrating  fiber  breakage  with  matrix  cracking  and 
degradation  are  shown  in  Figure  14.  The  calculation  was  preformed  with  the  post¬ 
failure  degradation  rules  implemented  in  DYNA3D,  i.e.  instantaneous  unloading. 
This  unloading  behavior  Weis  previously  shown  in  Figures  11  and  12.  The  calcu¬ 
lation  wa.s  preformed  with  a  three- layer  cross-ply  laminate,  [90°,  0°,  90°],  where  0° 
corresponds  to  the  x  direction  and  90°  corresponds  to  the  y  direction  of  the  laminate, 
as  shown  in  Figure  15.  A  monotonically  increasing  strain  was  applied  to  the  laminate 
in  the  x  direction,  i.e.  the  tensile  strain  was  applied  in  the  fiber  direction  of  the  0° 
lamina  and  in  the  transverse  fiber  direction  of  the  90°  laminae.  A  monotonically 
increasing  shear  strain  was  applied  simultaneously.  The  laminate  was  2.54-mm-thick 
with  the  following  tensile  strengths: 


Xr  = 

1379 

MPa 

(200  ksi) 

Yt  = 

34 

MPa 

(  5  ksi) 

Yc  = 

69 

Mpa 

(  10  ksi) 

S  = 

83 

MPa 

(  12  ksi) 

‘The  constant-stress  post-failure  degradation  rules  implemented  in  DYNA3D  were  modified  by 
APTEK  to  allow  the  shear  and/or  transverse  stresses  to  degrad  instantaneously  following  matri.x 
failure. 


(a)  Laminate  stress  resultant-strain  behavior  in  x  direction. 


(b)  Laminate  shear  stress  resultant-strain  behavior. 
Figure  14.  Demonstration  of  post-failure  behavior  in  DYNA.3D. 
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Figure  15.  Cross-ply  laminate  used  in  post-failure  demonstration  calculations. 


The  first  failure  to  occur  is  matrix  cracking  in  both  90°  laminae.  Once  these  lamina 
fail,  the  effective  modulus  of  the  laminate  is  reduced  because  the  failed  laminae  cannot 
carry  any  stress  in  the  transverse  fiber  direction.  The  change  in  effective  modulus 
of  the  laminate  is  indicated  by  the  drop  in  the  in-plane  stress  resultant  in  the  x 
direction,  N^,  when  Cj.  =  0.46%.  The  effective  shear  modulus  of  the  laminate  is 
also  reduced  to  one-third  its  original  value  because  the  failed  lamina  cannot  carry 
any  shear  stress  once  matrix  cracking  failure  occurs.  The  change  in  effective  shear 
modulus  is  indicated  by  the  drop  in  the  shear  stress  resultant,  when  e^y  —  0.46%. 

After  failure  of  the  90°  laminae,  stress  in  the  x  direction  and  shear  stress  must 
be  carried  by  the  0°  lamina.  However,  all  lamina  are  still  able  to  carry  stress  in 
the  y  direction.  As  the  laminate  extends  in  the  x  direction,  it  contracts  in  the  y 
direction  due  to  the  effective  Poisson’s  ratio  of  the  laminate.  This  ratio  depends  on 
the  elastic  constants  and  orientations  of  all  plies  within  the  laminate  and  is  calculated 
from  the  extensional  stiffness  coefficients  of  the  laminate;  i/^y  =  Ai2/A22^-  For  the 
cross-ply  laminate  used  in  this  example,  the  effective  Poisson’s  ratio  of  the  laminate. 
Uj;y  =  0.019,  is  less  than  the  Poisson’s  ratio  of  the  0°  lamina,  i/i2  =  0.3.  This  mismatch 
between  the  Poisson’s  ratios  of  the  laminate  and  the  0°  lamina  causes  tensile  stresses 
to  develop  in  the  transverse  fiber  direction  of  the  0°  lamina  and  compressive  stresses 
to  develop  in  the  fiber  direction  of  the  90°  lamina;  the  net  laminate  stress  resultant 
in  the  y  direction  is  zero.  .Matrix  cracking  failure  occurs  in  the  0°  lamina  when  the 
transverse  tensile  stress  exceeds  the  corresponding  strength  of  the  material.  .At  this 
strain  level,  Cj.  =  e^-y  %  1.0%,  none  of  the  lamina  are  able  to  carry  shear  stress  so 
the  shear  stress  drops  to  zero.  The  effective  Poisson's  ratio  of  the  laminate  is  also 
reduced  to  zero. 

"Extensional  stiffness  coefficients  and  laminate  effective  moduli  are  discussed  in  a  previous  report 
to  DNA  in  Referenco  [28] 
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Figure  16.  Post-failure  degradation  behavior  in  fiber-composite  lamina. 


The  last  failure  to  occur  in  the  laminate  is  fiber  breakage  in  the  0°  lamina  B leak¬ 
age  occurs  when  the  stress  in  the  fiber  direction  of  the  0“  lamina  exceeds  the  longi¬ 
tudinal  tensile  strength  of  the  material. 

This  calculation  demonstrates  the  implementation  of  the  post-failure  degradation 
rules  in  DYNA3D.  We  are  not  sure  that  instantaneous  unloading  for  the  matrix  failure 
modes  is  appropriate  for  all  fiber-composite  materials.  Other  possible  methods  are 
constant  stress  after  matrix  failure  (transverse  and  shear)  or  gradual  decay  of  the 
transverse  stress  and  shear  stress  after  the  matrix  fails,  as  shown  schematically  is 
Figure  16  and  discussed  in  References  [17]  and  [19].  Another  method  (not  shown)  is 
partial  reduction  of  some  of  the  elastic  moduli,  as  discussed  in  Reference  [29]. 

2.2.2  Effect  of  Post-Failure  Behavior  on  Ultimate 
Strength. 

In  this  section  we  demonstrate  how  the  post-failure  behavior  of  the  matrix  affects 
the  ultimate  strength  and  stress  resultant-strain  history  of  a  layered  composite.  This 
is  done  by  comparing  the  calculated  response  of  a  cross-ply  laminate  exhibiting  con¬ 
stant  stress  after  matrix  failure  with  the  calculated  response  of  a  laminate  exhibiting 
instantaneous  unloading. 

The  response  of  a  three- layer  cross-ply  laminate  with  instantaneous  unloading 
of  the  transverse  and  shear  stresses  was  discussed  in  the  previous  subsection.  Now 
consider  the  same  laminate  and  applied  strain  history,  but  allow  the  transverse  and 
shear  stresses  to  remain  constant  after  matrix  cracking  failure.  Comparisons  of  stress 
resultant-strain  histories  calculated  for  the  two  post-failure  behaviors  are  shown  in 
Figure  17. 

As  previously  demonstrated,  the  first  failure  to  occur  is  matrix  cracking  in  the  90° 
lamina,  which  reduces  the  effective  axial  and  shear  moduli  of  the  laminate.  With  the 
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(a)  Laminate  stress  resultant-strain  behavior  in  x  direction. 


Shear  Strain  (percent) 

(b)  Laminate  shear  stress  resultant-strain  behavior. 


Figure  17.  Comparison  of  stress  resultant-strain  behaviors  for  laminates  with  different 
post-failure  behaviors. 
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constant  stress  degradation  rule,  the  changes  in  effective  moduli  of  the  laminate  are 
indicated  by  the  chajige  in  slope  of  the  corresponding  stress  resultant-strain  curves. 
The  matrix  continues  to  carry  the  transverse  and  shear  stresses  at  failure,  but  all 
additional  stresses  must  be  carried  by  the  0°  lamina. 

Additional  strain  was  applied  to  the  laminate  until  matrix  cracking  occured  in  the 
0°  lamina  due  to  the  mismatch  between  the  Poisson’s  ratios  of  the  ply  and  laminate. 
At  this  strain,  —  e.xy  ^  1.0%,  none  of  the  lamina  are  able  to  carry  additional  shear 
load  so  the  shear  stress  resultant  remains  constant  with  increasing  shear  strain  until 
fiber  breakage  occurs.  Fiber  breakage  occurs  when  a  combination  of  the  increasing 
longitudinal  stress  and  constant  shear  stress  exceeds  the  strength  of  the  0°  lamina  at 
=  ^xy  =  1.46%. 

Note  that  the  maximum  in-plane  stress  resultant  calculated  with  the  constant 
stress  unloading  assumption  is  only  Nx  =  699  kN/m  (3992  Ib/in).  This  value  is  40% 
lower  than  the  in-plane  stress  resultant  calculated  with  instantaneous  unloading  of 
the  matrix,  which  is  Nx  =  1169  kN/m  (6675  Ib/in).  The  constant  stress  unloading 
value  is  lower  because  the  matrix  continued  to  carry  shear  stress  after  matrix  failure, 
and  shear  stress  is  assumed  to  weaken  the  composite  and  contribute  to  fiber  breakage 
in  the  Chang  fiber  breakage  failure  criterion.  These  comparative  examples  indicate 
that  the  calculated  response  and  ultimate  strength  of  a  laminate  with  matrix  failure 
are  very  sensitive  to  the  post-failure  behavior  of  the  matrix. 

2.3  COMPARISON  OF  FAILURE  CRITERIA  WITH 
EXPERIMENTAL  RESULTS. 

In  this  section,  numerical  results  are  comp<ired  with  experimental  results  available  in 
the  literature  to  establish  the  validity  of  the  Chang  failure  criteria  implemented  in 
DYNA3D.  Three  sets  of  material  property  tests  are  analyzed:  lamina  off-axis  tests  to 
verify  the  tensile  failure  criteria  (fiber  and  matrix),  lamina  combined  stress  tests  to 
verify  the  matrix  failure  criteria,  and  biaxial  stress  tests  on  quasi- isotropic  laminates 
to  verify  the  fiber  failure  criteria. 

Results  presented  in  this  section  indicate  that  the  Chang  matrix  tensile  failure  cri¬ 
terion  is  in  fair  agreement  with  lamina  off-axis  test  data.  The  Chang  fiber  breakage 
criterion  is  also  in  fair  agreement  with  failure  strengths  measured  from  two  off-axis 
tests,  although  two  data  points  are  not  enough  to  validate  the  criterion.  Comparison 
of  the  fiber  failure  criterion  with  biaxial  stress  tests  on  quasi-isotropic  laminates  was 
inconclusive  in  validating  the  Chang  fiber  breakage  failure  criterion.  Comparison  of 
the  matrix  failure  criteria  with  combined  stress  test  data  indicate  that  the  Chang 
matrix  tensile  and  compressive  failure  criteria  implemented  in  DYNA3D  do  not  ade¬ 
quately  model  matrix  strength  under  combined  loading  conditions.  A  modified  form 
of  the  Tsai-Wu  criterion  is  in  better  agreement  with  the  test  data  and  is  therefore 
recommended  for  implementation  in  DYNA3D. 
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Figure  18.  Geometry  of  the  off-axis  tensile  specimen. 

2.3.1  Lamina  OfF-Axis  Tests. 

Lamina  tensile  strength  is  often  characterized  by  performing  uniaxial  tests  on  uni¬ 
directional  specimens  with  fibers  oriented  at-angle  to  the  loading  axis,  as  shown  in 
Figure  18.  The  measured  strength  depends  on  the  off-axis  angle,  9,  between  the  fiber 
direction  and  loading  axis.  Extensional-shear  coupling  occurs  when  fibers  are  oriented 
at-angle  to  the  loading  axis. 

The  theoretical  tensile  strength  of  the  specimen  is  determined  from  the  Chang 
fiber  breakage  and  matrix  cracking  criteria.  This  is  done  by  transforming  the  applied 
axial  stress,  into  stresses  in  the  principal  material  directions  of  the  lamina,  as 
follows: 

(Ti 

0-2  =  (Tr  (15) 

(Ti2  —mn 

where  m  =  cos  6,  n  =  sin0,  and  9  is  the  orientation  angle  between  the  fiber  direction 
and  loading  axis.  The  state  of  stress  in  the  fiber  coordinate  system  is  biaxial.  The 
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maximum  applied  stress  is  found  by  substituting  the  lamina  stresses  from  Equation  15 
into  the  fiber  breakage  and  matrix  cracking  criteria,  Equations  12  and  13  respectively: 


Fiber  Breakage  :  <Tj.  Ifita-  = 


Matrix  Cracking  :  (Tx  jnuurix  = 


(16) 

(17) 


Equations  16  and  17  give  the  ultimate  tensile  strengths  for  fiber  breakage,  Cr  Ifiber- 
and  matrix  cracking,  Imatrixt  versus  off-axis  angle.  These  equations  are  plotted  in 
Figure  19  for  a  graphite-epoxy  specimen. 


For  small  angles,  the  ultimate  strength  calculated  for  fiber  breakage  is  less  than 
that  for  matrix  cracking.  For  larger  angles,  the  ultimate  strength  calculated  for  matrix 
cracking  is  less  than  that  for  fiber  breakage.  The  orientation  angle  which  separates 
the  fiber  and  matrix  modes  is  found  by  setting  the  ultimate  strength  for  fiber  breakage 
equal  to  the  ultimate  strength  for  matrix  cracking,  cTj.  lfiber=<^r  Imatnxi  which  gives  the 
following  separation  angle: 


(18) 


This  separation  angle  is  ^  =  8.8°  for  graphite-epoxy  with  Xt  =  1980  MPa  and 
Yt  =  48  MPa. 


Strengths  measured  from  an  off-axis  test  on  a  graphite-epoxy  specimen  are  com¬ 
pared  with  the  Chang  tensile  failure  criteria  in  Figure  19.  The  data  were  obtained 
from  Reference  [30,  p.  115].  There  is  good  agreement  between  the  test  data  and  the 
tensile  failure  criteria  implemented  in  DYNA3D.  The  fiber  breakage  criterion  is  in 
agreement  with  the  test  data  for  angles  less  than  the  separation  angle  while  the  ma¬ 
trix  cracking  criterion  is  in  agreement  with  the  test  data  for  angles  greater  than  the 
separation  angle.  Although  the  mode  of  failure  was  not  reported  for  the  test  data, 
the  comparisons  suggest  that  the  tests  were  terminated  once  matrix  tensile  failure 
occured  in  the  lamina  for  angles  greater  than  the  separation  angle. 

The  test  data  also  indicate  that  the  ultimate  strength  of  the  off-axis  specimen 
tested  at  5°  is  60%  less  than  a  specimen  tested  at  0°.  However,  this  test  does  not 
validate  the  fiber  breakage  failure  criterion  in  DYNA3D  for  two  reasons.  First,  the 
observed  failure  mode  was  not  reported  for  the  data,  so  we  do  not  know  if  fiber  or 
matrix  failure  occured  at  5°.  Second,  longitudinal  fiber,  transverse  fiber,  and  shear 
stresses  act  on  the  composite  in  the  fiber  coordinate  system,  but  only  shear  and 
longitudinal  stress  are  assumed  to  contribute  to  fiber  breakage  in  the  Chang  failure 
criterion.  Shear  stress  is  assumed  to  cause  fiber-matrix  debonding,  thus  reducing 
the  longitudinal  tensile  strength.  However,  stress  in  the  transverse  fiber  direction 
could  also  cause  matrix  cracking  failure,  which  may  also  contribute  to  a  reduction  in 
longitudinal  tensile  strength.  The  relative  effects  of  the  transverse  (crji)  and  shear 
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Figure  19.  Ultimate  tensile  strength  versus  ofF-axis  angle  for  unidirecti<^na'  tensile  test 
specimens. 

((7i2)  stresses  on  the  reduction  in  longitudinal  tensile  strength  are  not  established. 
Note:  it  could  also  be  argued  that  the  matrix  cracking  criterion  is  not  adequately 
verified  because  the  relative  effects  of  <tu  and  (T^  on  matrix  cracking  failure  are  not 
established. 

The  same  graphite-epoxy  test  data  are  compeired  with  two  commonly  used  failure 
criteria,  the  Tsai-Wu  and  maximum  stress  criteria,  in  Figure  20.  This  was  done  to 
determine  how  the  form  of  the  failure  criteria  affects  the  shape  of  ultimate  strength 
versus  off-axis  angle  curve.  The  Tsai-Wu  criterion  is  a  single  smooth  curve  in  ultimate 
stress  versus  orientation-angle  space,  and  is  in  excellent  agreement  with  the  data. 
The  maximum  stress  criteria  consists  of  three  intersecting  curves  which  correspond 
to  longitudinal  fiber  tensile  failure,  transverse  tensile  failure,  and  shear  failure.  The 
maximum  stress  criteria  are  also  in  good  agreement  with  the  data,  although  the 
maximum  stress  criteria  overestimate  the  strength  of  the  laminate  at  .30°  because 
interaction  between  the  in-plane  shear  and  transverse  stress  components  are  not  taken 
into  account.  In  conclusion,  differences  between  the  Chang,  Tsai-Wu.  and  maximum 
stress  criteria  are  not  readily  apparent  when  the  failure  curves  are  plotted  in  ultimate 
strength  versus  orientation-angle  space. 

Another  way  to  compare  off-axis  test  data  with  failure  criteria  is  in  lamina  stress 
(  <^11  1  <^22  1  <^12)  space.  This  comparison  provides  a  more  exacting  assessment  of  the 
criteria  than  the  ultimate  strength  vs.  orientation  angle  plots,  particularly  for  inter¬ 
mediate  orientation  angles  where  failure  occurs  in  the  matrix  and  is  dominated  by 
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Orientation  Angle  (degrees) 

Figure  20.  Comparison  of  the  Tsai-Wu  and  maximum  stress  criteria  with  data  from 
off-axis  test  specimens. 

<Ti2  and/or  <722-  These  stresses,  and  cr22,  are  small  compared  to  the  maximum 
ordinate  of  the  ultimate  stress  vs.  orientation  angle  plot,  which  is  cr^  ss  Xt- 

For  the  Chang  criteria,  comparisons  in  stress  space  are  easy  to  visualize  because 
failure  by  fiber  breakage  or  matrix  cracking  are  each  associated  with  surfaces  parallel 
to  one  of  the  stress  axes  (cylindrical  surfaces),  defined  by  Equations  12  and  13.  The 
only  difficulty,  as  previously  mentioned,  is  that  one  must  know  which  mode  each 
specimen  failed  in  to  plot  its  data  point  for  comparison  with  the  appropriate  failure 
surface.  Comparisons  are  shown  in  Figures  21  and  22  with  the  assumption  that 
failure  occured  by  fiber  breakage  for  0  <  8.8°,  and  by  matrix  cracking  for  0  >  8.8°. 
Figures  21  and  22  can  be  visualized  looking  down  the  axis  of  each  cylindrical  failure 
surface.  Thus,  the  distance  from  a  test  data  point  to  the  surface  is  seen  as  a  true 
length.  To  help  interpret  these  figures,  measured  failure  stresses  in  the  principal 
material  directions  of  the  lamina  are  listed  in  Table  4. 

Figure  21  shows  fiber  failure  data  plotted  in  a  ,  cr]2  projection,  with  the  associ¬ 
ated  end  view  of  tiie  Chang  fiber  breakage  cylindrical  failure  surface.  The  data  point 
for  (712  =  0  (^  =  0)  lies  exactly  on  the  surface  because  the  data  from  this  test  was 
used  in  establishing  the  surface.  The  shear  strength  of  85  MPa  at  dn  =  0  completed 
the  surface  definition  and  was  obtained  from  tensile  test  data  on  ±45°  lamina  spec¬ 
imens  (Reference  [30],  p.  63).  The  other  fiber  failure  point  from  the  off-axis  tests  at 
0  =  .5°  lies  below  the  surface  by  a  more  readily  apparent  distance  than  in  the  ultimate 
strength  plot  in  Figure  19. 
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Table  4.  Graphite-epoxy  off-axis  test  data. 
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Figure  21.  Comparison  of  off-axis  test  data  with  Chang  fiber  breakage  criterion. 
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Figure  22.  Comparison  of  off-axis  test  data  with  Chang  matrix  cracking  criterion. 

Figure  22  shows  matrix  failure  data  plotted  in  a  <722 ,  <712  projection  along  with 
the  associated  end  view  of  the  Chang  matrix  failure  cylindrical  surface.  Again,  the 
data  point  for  <712  =  0  lies  exactly  on  the  surface,  because  it  was  used  to  establish 
the  surface.  The  important  observation  here  is  that  the  data  point  measured  at  30° 
(<722  =  28  MPa,  <7i2  =  48  MPa)  falls  substantially  off  the  failure  surface  when  viewed 
along  the  ctu  stress  axis,  because  the  large  value  of  <73,  =  Xt  does  not  dominate  the 
comparison  as  it  does  in  the  ultimate  strength  plot  in  Figure  19. 

Conclusions. 

The  Chang  fiber  breakage  and  matrix  cracking  criteria  are  only  in  fair  agreement 
with  the  test  data  when  the  data  are  plotted  in  stress  space.  The  two  data  points 
in  Figure  21  are  not  enough  to  validate  the  fiber  breakage  criterion.  In  addition, 
the  effects  of  <722  on  fiber  breakage  failure  and  (Tu  on  matrix  cracking  failure  are 
not  adequately  addressed  in  these  tests.  Additional  biaxial  characterization  data  of 
composite  specimens  is  needed  to  validate  both  the  Chang  fiber  and  matrix  tensile 
failure  criteria.  Additional  biaxial  tests  are  discussed  in  the  following  two  subsections. 

2.3.2  Fiber  Failure  Modes  Under  Biaxial  Stress. 

Swanson  et  al.  in  References  [3 1 -.34]  reported  biaxial  stress  experiments  on  quasi¬ 
isotropic  carbon-epoxy  (AS4/3501-6)  laminates,  [90°,  ±45°,  0°]s.  Fiber  failure  data 
were  obtained  from  static  tests  on  thin-walled  cylindrical  tubes  under  combined  in¬ 
ternal  pressure  and  axial  load  (both  tension  and  compression).  The  failure  data  are 
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Figure  23.  Failure  strain  data  for  a  quasi-isotropic  carbon-epoxy  laminate  under  bi¬ 
axial  stress. 

shown  in  Figure  23  for  combinations  of  hoop  and  axial  strain.  The  average  hoop  fail¬ 
ure  strain  from  all  tubes  tested  in  tension-tension  (axial  tension  combined  with  hoop 
tension)  is  1.40%.  This  average  failure  strain  is  close  to  a  value  of  1.46%  measured 
by  Swanson  in  unidirectional  tensile  coupon  tests.  Swanson  et  al.  concluded  that  the 
tensile  hoop  failure  strain  is  constant  at  failure  and  independent  of  the  applied  stress 
state,  at  lea^t  in  the  tension- tension  stress  quadrant. 

Swanson  et  al.  in  Reference  [32]  also  reported  combinations  of  hoop  and  axial 
stress  which  produced  failure  in  the  laminate.  They  compared  the  failure  stress  data 
with  failure  stress  envelopes  calculated  using  the  maximum  fiber  strain  criterion  and 
the  Tsai-Wu  failure  criterion.  These  failure  envelopes  are  compared  with  the  test 
data  in  Figure  24  for  the  tension-tension  stress  quadrant  only.  Two  envelopes  are 
shown  for  the  maximum  strain  criterion.  The  dcished  envelope  is  calculated  with 
linear  laminated  plate  theory  (LPT)  and  the  solid  v-shape  line  is  calculated  with  a 
progressive  softening  model  that  includes  nonlinear  shear  behavior  and  reduction  in 
ply  stiffness^  due  to  matrix  cracking.  The  maximum  strain  criterion  with  progressive 
softening  agrees  best  with  the  measured  data,  although  the  maximum  strain  criterion 
without  the  degradation  model  is  also  in  re^lsonable  agreement  with  the  data.  The 
Tsai-Wu  criterion  is  in  poor  agreement  with  the  data.  With  the  Tsai-Wu  criterion. 

^Swanson  and  Christoforou  in  Reference  [32]  use  a  gradual  degradation  rule  to  model  the  reduc¬ 
tion  in  transverse  Young’s  modulus,  £’22,  as  a  function  of  transverse  strain. 
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Figure  24.  Laminate  failure  envelopes  calculated  with  the  Tsai-Wu  and  maximum 
strain  criteria  for  a  quasi-isotropic  carbon-epoxy  laminate  under  biaxial 
stress. 

‘last  ply’  laminate  failure  occurs  when  the  Tsai-Wu  criterion  is  satisfied  for  all  plies. 
These  results  suggest  that  fiber  failure  is  not  adequately  modeled  by  the  Tsai-Wu 
criterion  under  biaxial  stress  conditions. 

Swanson  et  al.  did  not  report  whether  fiber  failure  initiated  in  the  hoop,  axial, 
or  at-angle  plies.  Therefore,  the  test  data  are  interpreted,  in  the  present  report, 
by  calculating  the  failure  envelopes  for  the  individual  lamina.  Failure  envelopes  are 
calculated  for  two  different  failure  criteria:  the  maximum  stress  criterion,  and  an 
interactive  failure  criterion.  The  failure  envelope  for  the  entire  laminate  is  determined 
from  the  failure  envelopes  for  the  individual  lamina.  Results  from  these  caJculations 
are  discussed  in  the  following  paragraphs. 

Maximum  Stress  Failure  Criteria. 

Combinations  of  hoop  and  axial  stress  which  satisfy  the  maximum  stress  criterion 
for  fiber  failure  are  shown  in  Figure  25.  For  each  lamina,  fiber  failure  is  bounded 
by  parallel  lines  in  stress  space.  For  example,  two  straight  lines  are  shown  for  the 
axial  lamina:  one  for  fiber  tensile  failure  and  the  other  for  fiber  compression  failure. 
The  semi-infinite  region  bounded  by  these  lines  represents  the  stress  states  the  axial 
lamina  can  sustain  without  fiber  tensile  or  compressive  failure  occuring.  The  hoop 
and  ±45°  laminae  fail  in  fiber  compression  outside  the  range  of  this  figure. 

Failure  in  the  laminate  is  the  minimum  interaction  diagram  of  the  individual  plies. 
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Figure  25.  Laminate  failure  envelope  calculated  with  the  maximum  stress  criterion 
for  a  quasi-isotropic  carbon-epoxy  laminate  under  biaxial  stress. 


The  shaded  area  represents  the  biaxial  stress  states  which  the  laminate  can  sustain 
without  fiber  failure  occurring  in  any  lamina.  The  shaded  area  is  bounded  by  the  fiber 
tensile  failure  curve  for  the  hoop  lamina,  and  the  fiber  tensile  and  compressive  failure 
curves  for  the  axial  lamina.  The  laminate  failure  envelope  is  in  good  agreement  with 
the  test  data,  shown  as  dots  in  Figure  25,  reported  by  Swanson  et  al. 

Note:  The  following  lamina  properties  and  strengths  were  reported  by  Swanson  in 
Reference  [32]  and  used  in  the  calculations:  En  =  131  GPa,  E22  =  11.2  GPa,  1/12  — 
0.28,  Xt  =  1950  MPa.  Lamina  compressive  strength  and  in-plane  shear  modulus 
and  strength  were  not  reported.  A  compressive  strength  of  Xc  =  1074  MPa  was 
estimated  from  the  measured  biaxial  stress  data.  A  shear  modulus  of  G12  =  6600 
MPa  was  calculated  using  laminated  plate  theory,  the  lamina  elastic  constants,  and 
laminate  stiffness  measurements  reported  in  Reference  [32].  A  shear  strength  of  5  = 
52  MPa  was  estimated  from  data  for  AS4/55A  carbon-epoxy,  which  was  reported  in 
Reference  [35].  Reference  [36]  reports  shear  strengths  ranging  from  about  40  MPa  for 
T300/5208  carbon-epoxy  tape  to  125  MPa  for  AS4/3502  carbon-epoxy  tape  at  room 
temperature. 

Interactive  Failure  Criterion. 

Failure  envelopes  were  also  calculated  using  an  interactive  criterion  of  the  following 
form: 
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where  X  is  the  longitudinal  fiber  strength  in  either  tension  or  compression.  For  fiber 
tensile  failure  (an  >  0),  X  =  Xj,  and  Equation  19  reduces  to  the  Chang  fiber  break¬ 
age  criterion.  For  fiber  compressive  failure  (an  <  0),  X  =  Xc  in  Equation  19;  this 
fiber  compressive  failure  criterion  is  suggested  by  Humphreys  and  Rosen  in  Refer¬ 
ence  [37]. 

Combinations  of  hoop  and  axial  stress  which  satisfy  the  interactive  criterion  for 
fiber  failure  are  shown  in  Figure  26.  For  the  hoop  and  axial  lamina,  fiber  failure  is 
bounded  by  parallel  lines  in  stress  space,  in  agreement  with  those  calculated  from  the 
maximum  stress  criterion.  This  is  because  the  shear  stress  in  the  principal  material 
directions  of  the  hoop  and  axial  lamina  is  zero.  Fiber  failure  occurs  when  the  stress  in 
the  longitudinal  fiber  direction  exceeds  the  corresponding  strength,  i.e.  the  interactive 
criterion  reduces  to  the  maximum  stress  criterion  when  the  shear  stress  is  zero. 

For  the  ±45°  plies,  the  shear  stress  in  the  principal  material  directions  is  not 
zero.  Hence  the  lamina  failure  envelope  calculated  with  the  interactive  criterion. 
Equation  19,  is  an  ellipse  in  stress  sp  ice.  Tensile  or  compressive  fiber  failure  occurs 
in  the  at-angle  plies  for  stress  states  outside  the  region  bounded  by  the  ellipse.  Note 
that  the  elliptical  failure  envelope  for  the  at-angle  plies  is  contained  almost  entirely 
within  the  region  bounded  by  the  hoop  and  axial  lamina.  The  portion  of  the  elliptical 
failure  envelope  that  extends  outside  the  bounds  for  tensile  failure  in  the  hoop  and 
axial  lamina  is  shown  as  a  dashed  line. 

Failure  in  the  laminate  is  the  minimum  interaction  diagram  of  the  individual  plies. 
The  shaded  area  is  bounded  by  the  laminate  failure  envelope.  Outside  the  shaded 
area  are  the  stress  states  for  which  fiber  tensile  or  compressive  failure  occured  in 
at  least  one  lamina.  Note  that  this  area  is  bounded  almost  entirely  by  the  failure 
envelopes  for  the  ±45°  plies  (the  ellipse). 

The  laminate  failure  envelope  calculated  from  the  interactive  failure  criterion  is  not 
in  good  agreement  with  the  test  data.  The  calculated  failure  strength  of  the  laminate 
is  strongly  influenced  by  the  shear  stress  and  strength  of  the  at-angle  lamina.  It 
is  possible  that  the  shear  strength  used  in  the  calculation  is  not  in  agreement  with 
that  of  the  test  specimens,  e.g.  a  shear  strength  of  100  MPa  (not  shown)  gives  better 
agreement  with  the  data.  It  is  also  possible  that  matrix  cracking  occured  in  the  at- 
angle  plies,  which  would  reduce  the  value  of  the  shear  modulus,  and  hence  reduces 
the  tendency  for  fiber  breakage  in  the  ±45°  plies.  Swanson  et  al.  reported  that  the 
measured  stress-strain  response  of  the  laminate  shows  a  small  reduction  in  laminate 
stiffness  at  higher  stress  levels,  which  they  attributed  to  matrix  failure.  The  results 
from  a  sample  calculation  with  the  shear  modulus  reduced  to  half  its  original  value 
are  shown  in  Figure  27.  With  Gu  =  3300  MPa,  the  laminate  failure  envelope  is  in 
better  agreement  with  the  test  data.  With  the  shear  modulus  reduced  to  Gu  =  0. 
the  failure  envelope  calculated  with  the  maximum  stress  failure  criterion  is  obtained. 


40 


Figure  26.  Laminate  failure  envelope  calculated  with  an  interactive  failure  criterion 
for  a  quasi-isotropic  carbon-epoxy  laminate  under  biaxial  stress. 


Figure  27.  Laminate  failure  envelope  calculated  with  an  interactive  failure  criterion 
and  a  reduced  shear  modulus. 
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Conclusions. 


Comparisons  between  the  test  data  and  laminate  failure  envelopes  in  Figure  26  and 
27  are  not  definitive  in  validating  the  interactive  failure  criterion:  the  comparisons 
neither  refute  nor  support  the  data.  The  comparisons  do  not  refute  the  interactive 
criterion  if  the  shear  modulus  is  reduced  to  account  for  matrix  cracking  failure  in  the 
at-angle  plies.  On  the  other  hand,  the  comparisons  do  not  support  the  interactive 
criterion  because  the  contribution  to  fiber  failure  from  shear  stress  is  reduced  by 
matrix  cracking,  so  the  full  effect  of  the  shear  stress  term  in  Equation  19  is  not 
validated.  In  fact,  the  maximum  stress  criterion,  which  neglects  the  contribution  of 
shear  stress  to  fiber  failure,  is  also  in  good  agreement  with  the  test  data.  Additional 
test  data  are  needed  to  validate  the  interactive  criterion,  such  as  biaxial  test  data  on 
unidirectional  specimens  under  combined  longitudinal  fiber  stress  and  shear  stress. 
In  particular,  we  want  to  validate  the  Chang  fiber  breakage  criterion'*  implemented  in 
DYNA3D,  and  validate  a  fiber  compressive  failure  criterion  for  future  implementation 
into  DYNA3D. 

The  analyses  completed  to  date  also  indicate  the  need  to  accurately  establish 
the  reduction  in  lamina  properties  with  matrix  cracking  failure,  particularly  the 
post-failure  implementation  in  DYNA3D.  Comparisons  of  the  DYNA3D  post-failure 
degradation  rules  with  experimental  results  available  in  the  literature  are  discussed 
in  Section  2.4. 

2.3.3  Matrix  Failure  Modes  Under  Combined  Stress. 

Swanson,  Messick,  and  Tian  in  Reference  (35j  reported  failure  strengths  of  carbon- 
epoxy  lamina  under  conditions  of  combined  transverse-fiber  and  shear  stress.  Com¬ 
bined  torsion  and  axial  load  (both  tension  and  compression)  tests  were  conducted 
using  unidirectional  hoop- wound  tubes.  Tube  tests  help  to  minimize  free  edge  de¬ 
lamination  and  end  condition  (gripping  method)  effects  which  are  present  in  flat  test 
specimens. 

Matrix  failure  data  for  carbon-epoxy  are  shown  in  Figure  28.  Transverse  tensile 
and  compressive  failure  stresses  are  plotted  as  a  function  of  shear  stress.  The  data 
show  a  strong  interaction  between  strength  in  the  transverse  fiber  direction  and  shear 
strength.  Note  that  the  shear  strength  measured  with  moderate  transverse  compres¬ 
sive  stress  is  greater  than  the  pure  shear  strength:  the  shear  strength  is  approximately 
73  MPa  for  <732  =  —50  MPa  while  the  shear  strength  is  52  MPa  for  (T22  =  0  MPa. 

Combinations  of  transverse  fiber  and  shear  stress  which  satisfy  the  Chang  matrix 

■*  Analysts  who  use  DYNA3D  should  be  aware  of  the  extent  to  which  shear  stress  contributes  to 
fiber  breakage  in  the  Chang  criterion.  Currently,  fiber  breakage  is  calculated  in  DYNA3D  when 
a  pure  shear  load  is  applied  to  a  laminated  shell  element.  However,  the  criterion  is  intended  to 
calculate  fiber  breakage  in  composites  loaded  primarily  in  the  longitudinal-fiber  direction.  Shear 
stress  contributes  to  fiber  breakage  by  weakening  the  fiber-matrix  interface  or  through  shear  failure 
in  the  matrix. 
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Figure  28.  Comparison  of  Chang  matrix  failu-e  criteria  with  matrix  failure  data  for 
carbon-epoxy. 

tensile  and  compressive  failure  criteria®  are  also  plotted  in  Figure  28.  The  calculated 
tensile  failure  curve  overestimates  the  tensile  strength  under  combined  load  condi¬ 
tions.  The  calculated  compressive  failure  curve  is  in  fair  agreement  with  the  test  data 
for  laxge  transverse  compressive  stresses  but  does  not  predict  the  increase  in  shear 
strength  with  moderate  compressive  stress  observed  in  the  test  data. 

Because  of  the  poor  agreement  between  calculated  and  measured  strengths  in  Fig¬ 
ure  28,  a  study  was  performed  to  determine  if  other  combinations  of  Yc  and  S  provide 
better  agreement.  The  Chang  matrix  compressive  failure  criterion®.  Equation  14.  is 
repeated  here  to  aid  the  discussion: 


The  matrix  compression  failure  criterion  contains  a  linear  stress  term  whose  sign 
depends  on  the  magnitude  of  the  uniaxial  compressive  strength,  Yc,  and  the  shear 
strength,  5.  Equation  20  describes  an  ellipse  in  stress  space:  one-fourth  of  the  ellipse 
is  shown  in  Figure  29  for  three  combinations  of  Yc  and  5.  When  Yc  =  25,  the 
coefficient  of  the  linear  stress  term  is  zero  and  the  axes  of  the  ellipse  coincide  with  the 
coordinate  axes  (<722  and  cri2  axes).  The  maximum  shear  stress  occurs  when  <722  =  0. 

^Uniaxial  strengths  and  properties  from  Reference  [35]  were  used  to  calculate  the  failure  surfaces: 
Yt  =  27  Mpa,  Yc  =  95  MPa,  and  S=52  MPa,  G12  =  4300  MPa,  and  q  =  6.88e  -  08  .MPa^. 

^The  nonlinear  shear  stress  parameter,  a,  has  been  set  to  zero  in  Equation  14  to  simplify  the 
analysis. 
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Figure  29.  Chang  matrix  compression  failure  envelope  calculated  from  three  combina¬ 
tions  of  transverse  compressive  strength,  Vc,  in-plane  shear  strength, 

5. 

When  Vc  7^  25,  one  axis  of  the  ellipse  is  parallel  to  the  (T12  coordinate  axis,  but  does 
not  coincide  with  the  cri2  axis.  When  Yc  <  25,  the  maximum  shear  stress  occurs 
when  the  transverse  stress  is  tensile  (not  shown).  When  Yc  >  25,  the  maximum  shear 
stress  occurs  when  the  transverse  stress  is  compressive.  This  behavior,  i.e.  an  increase 
in  shear  strength  with  transverse  compressive  stress,  is  observed  in  the  measured  data 
in  Figure  28. 

The  maximum  shear  stress  is  determined  by  maximizing  cri2  in  Equation  20  with 
respect  to  a’22-  The  maximum  shear  strength  under  combined  loading,  5n,j»x>  is: 


9  -  9 

‘^rriAX  —  ^ 


1 

2^ 


(21) 


The  shear  strength  is  a  maximum  when  the  compressive  stress  is; 


<722 


[K^-4521 

2Yc 


(22) 


Equations  21  and  22  are  valid  for  any  combination  of  Yc  and  5.  Note  that  when 
Yc  »  25,  the  maximum  shear  stress  attains  its  largest  values  as  <722  =t'  —  Vc/2  . 

The  measured  shear  strength  under  combined  loading  in  Figure  28  is  about  1.4.') 
times  the  pure  shear  strength.  Substitution  of  5niax  =  1.45  5  into  Equation  21  gives 
Yc  —  5,9.  This  combination  of  Yc  and  5  (not  shown)  still  does  not  fit  the  measured 
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data  well.  Although  the  Chang  matrix  compression  failure  criterion  allows  for  in¬ 
creased  shear  strength  with  transverse  compressive  stress,  a  good  fit  to  the  measured 
data  cannot  be  obtained  with  any  reasonable  combinations  of  Yc  and  S.  Two  mod¬ 
ifications  to  the  matrix  compression  failure  criterion  implemented  in  DYNA3D  are 
suggested  in  the  following  paragraphs. 

Hashin  Matrix  Compression  Failure  Criterion. 

One  modification  of  the  Chang  matrix  compression  failure  criteria  is  to  include  the 
transverse  shear  strength.  As  discussed  in  Section  1.3,  the  Chang  matrix  compres¬ 
sive  failure  criterion  is  a  modification  of  the  Hashin  criterion.  Hashin  derived  his 
two-dimensional  criterion  in  terms  of  both  the  in-plane  shear  strength,  5,  and  the 
transverse  shear  strength,  S't,  as  follows: 

This  criterion  is  easily  modified  to  include  nonfinear  shear  behavior  following  the 
method  of  Chang,  as  follows: 

(^)  ■^[(^) 

Transverse  shear  strength  is  not  reported  for  the  carbon-epoxy  data  in  Refer¬ 
ence  [35].  If  we  arbitrarily  assume  a  value  of  St  =  20  MPa,  the  Hashin  matrix 
compression  failure  criterion  (with  nonlineaj  sheax  behavior)  predicts  a  substantial 
increase  in  shear  strength  with  moderate  compressive  stress,  as  shown  in  Figure  30. 
This  example  demonstrates  that  including  the  transverse  shear  strength  in  the  matrix 
compression  failure  criterion  has  a  significant  effect  on  the  strength  of  the  composite 
under  combined  loading  conditions  if  St  is  substantially  less  than  5.  However,  we 
do  not  know  if  the  transverse  shear  strength  assumed  in  this  example  is  appropriate. 
Rosen  and  Hashin  in  Reference  [8]  state  that  the  transverse  shear  strength  is  difficult 
to  measure.  A  reasonable  approximation  for  Sj  is  the  ultimate  shear  strength  for  the 
epoxy  matrix  alone,  but  the  matrix  shear  strength  was  not  reported  in  Reference  [35]. 

Modified  Tsai-Wu  Failure  Criterion. 

Swanson  and  Christoforou  in  Reference  [31]  use  an  interactive  failure  criterion  that  is 
a  modified  form  of  the  Tsai-Wu  criterion.  The  criterion  models  both  the  transverse 
tensile  and  compressive  behavior  of  carbon-epoxy.  The  quadratic  criterion  contains 
a  linear  term  whose  sign  depends  on  the  magnitude  of  the  transverse  tensile  and 
compressive  strengths,  Yt  and  Yc,  as  follows: 

F2<T22  -h  F22<T22  +  P'66<^\2  ^  1  (25) 
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Figure  30.  Comparison  of  Hashin  matrix  failure  criteria  with  measured  data  for 
carbon-epoxy. 


where 

_  _  yc-rr 

-  -rPT 

F,,  =  ^  (26) 

F  -  ± 

This  criterion  fits  the  data  very  well,  as  is  shown  in  Figure  31.  Note  that  the  criterion 
is  obtained  by  neglecting  all  terms  containing  <7u  in  Equation  2.  The  effect,  if  any. 
of  (Til  on  the  matrix  failure  strength  is  an  issue  Swanson  et  al.  in  Reference  [35]  plan 
to  address  in  future  research. 


Conclusions. 

The  Chang  matrix  criteria,  the  present  criteria  in  DYNA3D,  do  not  adequately  model 
matrix  failure  observed  in  carbon-epoxy  composites.  When  test  data  becomes  avail¬ 
able.  additional  studies  should  be  performed  on  fiberglass-epoxy  and  Kevlar-epoxy 
composites  of  interest.  Based  on  the  studies  completed  to  date,  the  modified  Tsai- 
Wu  criterion  is  recommended  for  modeling  matrix  tensile  and  compressive  failure 
in  fiber  composites.  Also,  the  modified  Tsai-Wu  criterion  identifies  the  tensile  and 
compressive  matrix  failure  modes,  continuing  in  the  spirit  of  Chang  to  identify  fail¬ 
ure  modes,  so  that  progressive  damage  can  be  analyzed  by  including  degradation  of 
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Figure  31.  Comparison  of  modified  Tsai-Wu  criterion  with  matrix  failure  data  for 
carbon-epoxy. 


lamina  properties. 


Table  5.  Failure  criteria  currently  implemented  in  DYNA3D. 


Failure  Mode 

Present  Implementation 

Recommended  Implementation  or  Action 

Fiber  Tension 

Fiber  Compression 

Matrix  Tension 

Matrix  Compression 

Chang  Fiber  Breakage 

None 

Chang  Matrix  Cracking 
Chang  Matrix  Compression 

Verify  Chang  Fiber  Breakage  Criterion 
Verify  Interactive  Criterion  and 

Implement  into  DYNA3D 

Implement  Modified  Tsai-Wu  Criterion 
Implement  Modified  Tsai-Wu  Criterion 

2.3.4  Conclusions. 

Table  5  compares  the  failure  criteria  currently  implemented  in  DYNA3D  with  the 
criteria  recommended  for  implementation.  Based  on  comparisons  of  various  failure 
criteria  with  test  data  available  in  the  literature,  the  following  conclusions  are  made; 

•  Comparisons  of  both  the  maximum  fiber  tensile  stress  criterion  and  the  Chang 
fiber  breakage  criterion  with  carbon-epoxy  data  from  biaxial  stress  tests  were 
not  conclusive  in  validating  either  criterion.  This  is  because  the  maximum 
stress  criterion,  and  the  Chang  fiber  breakage  criterion  used  in  conjunction 
with  a  degraded  shear  modulus  to  account  for  matrix  cracking,  were  both  in 
reasonable  agreement  with  the  test  data.  However,  the  Tsai-Wu  criterion  *was 
in  poor  agreement  with  the  test  data,  which  suggests  that  fiber  failure  is  not 
adequately  modeled  by  this  criterion.  Additional  compeirisons  with  test  data  are 
needed  to  validate  the  criteria,  in  particular  the  Chang  fiber  breakage  criterion 
which  is  currently  implemented  in  DYNA3D. 

•  Comparisons  of  bom  the  maximum  fiber  compressive  stress  criterion  and  an 
interactive  fiber  compressive  failure  criterion  with  data  from  biaxial  stress  tests 
were  not  conclusive  in  validating  either  criterion.  A  fiber  compressive  failure 
criterion  is  not  currently  implemented  in  DYNA3D.  Additional  comparisons 
with  test  data  are  needed  to  validate  a  fiber  compressive  failure  criterion  for 
future  implementation  into  DYNA3D. 

•  The  modified  Tsai-Wu  criterion  for  matrix  tensile  and  compressive  failure  is  rec¬ 
ommended  for  implementation  into  DYNA3D  to  replace  the  Chang  matrix  fail¬ 
ure  criteria.  The  modified  Tsai-Wu  criterion  is  in  good  agreement  with  failure 
strengths  of  carbon-epoxy  measured  under  conditions  of  combined  transverse- 
fiber  and  shear  stress.  The  Chang  matrix  failure  criteria  are  not  in  reason¬ 
able  agreement  with  the  test  data.  The  Chang  matrix  tensile  failure  criterion 
overestimates  the  matrix  tensile  strength  under  combined  load  conditions.  The 
Chang  matrix  compression  failure  criterion  does  not  predict  an  increase  in  shear 
strength  with  moderate  compressive  stress  comparable  with  that  observed  in  the 
test  data. 
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2.4  COMPARISON  OF  POST-FAILURE 

DEGRADATION  RULES  WITH  EXPERIMENTAL 
RESULTS. 

A  number  of  post-failure  material  property  degradation  rules  are  available  from  the 
literature,  see  for  example  the  survey  by  Nahas  in  Reference  [17].  Many  fit  into  one 
of  three  general  categories:  instantaneous  unloading,  gradual  unloading,  or  constant 
stress  after  ply  failure.  These  three  categories  were  previously  shown  in  Figure  14  of 
Section  2.2.  Behavior  of  the  ply  when  it  unloads,  as  well  as  which  elastic  constants  are 
degraded,  depends  on  the  mode  of  failure  of  the  particular  composite  being  analyzed. 

One  common  method  for  estimating  stiffness  reductions  in  laminates  with  cracked 
plies  is  the  ply  discount  method.  With  this  method,  one  or  more  of  the  elastic  constants 
of  a  cracked  ply  are  set  equal  to  zero.  This  method  belongs  to  the  instantaneous 
unloading  category.  For  example,  if  matrix  cracking  failure  occurs  in  a  ply  then  the 
following  are  set  to  zero:  Young’s  modulus  in  the  transverse-fiber  direction.  E22,  the 
in-plane  shear  modulus,  Gu,  and  the  major  and  minor  Poisson’s  ratio’s.  1/12  and  1/21. 
The  ply  is  still  able  to  carry  longitudinal  load  in  the  fiber  direction,  i.e.  ^  0.  It 
is  also  assumed  that  degradation  is  restricted  to  the  ply  that  cracks.  The  stiffness 
coefficients  of  the  laminate  are  recalculated  after  the  properties  of  the  cracked  ply  are 
reduced.  A  ply  discount  scheme  is  presently  implemented  in  DYNA3D. 

Numerous  studies  are  available  from  the  literature  for  determining  the  stress- strain 
response  of  laminates  with  ply  failure  due  to  matrix  cracking.  Some  of  the  earliest 
studies  were  done  by  Petit  and  Waddoups  in  Reference  [38],  who  suggest  a  gradual 
failure  model  for  boron-epoxy  laminates,  and  Hedm  and  Tsai  in  Reference  [39],  who 
suggest  a  constant  stress  model  for  graphite-epoxy  and  glass-epoxy  cross-ply  lami¬ 
nates.  More  recent  studies  on  transverse  cracking  and  laminate  stiffness  reductions 
were  performed  by  Talreja  in  Reference  [40]  and  Highsmith  and  Reifsnider  in  Ref¬ 
erence  [41].  These  studies  are  discussed  in  this  section.  Few  studies  are  available 
for  determining  laminate  stress-strain  response  in  compression.  Petit  and  Waddoups 
in  Reference  [38]  suggest  using  a  gradual  failure  model  for  ply  matrix  compression 
failure. 

We  begin  this  section  with  a  description  of  matrix  crack  formation.  This  is  fol¬ 
lowed  by  a  discussion  of  results  presented  by  Talreja,  and  Highsmith  and  Reifsnider. 
on  stiffness-reduction  mechanisms  in  cracked  laminates.  In  the  final  section,  modifi¬ 
cations  to  the  post-failure  degradation  rules  currently  implemented  in  DYNA3D  are 
recommended. 

2.4.1  Micromechanics  of  Matrix  Cracking. 

The  initiation  and  subsequent  multiplication  of  matrix  cracks  in  transverse  or  off- 
axis  plies  of  laminates  subjected  to  longitudinal  tension  is  presented  in  this  section. 
Cracks  initiate  in  transverse  or  off-axis  plies  at  relatively  low  values  of  applied  stress 
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(a)  STATE  0  (b)  STATE  1  (c)  STATE  2 

Figure  32.  Multiplication  of  matrix  cracks  in  the  transverse  ply  of  a  laminate  loaded 
in  longitudinal  tension. 

or  strain,  25%  to  33%  of  the  matrix  fracture  stress,  as  discussed  in  Reference  [42]. 
Cracks  initiate  at  areas  where  flaws  exists  or  where  fibers  debond  from  the  matrix. 
The  initial  cracks  are  widely  spaced. 

Nuismer  and  Tan,  in  Reference  [43],  developed  a  theoretical  model  in  which  crack 
multiplication  proceeds  in  a  deterministic  manner.  Consider  a  symmetric  cross-ply 
laminate  with  an  interior  90“  ply  (transverse  to  the  applied  strain)  containing  a  pe¬ 
riodic  array  of  cracks  of  spacing  2L,  as  shown  in  Figure  32a.  As  the  applied  strain  is 
slowing  increaised  from  ^  to  Ci,  a  new  strain  level  is  reached  at  which  a  new  set  of 
cracks  is  formed.  The  crack  formation  is  Msumed  to  be  deterministic  with  the  newly 
formed  cracks  equally  spaced  between  the  original  cracks.  This  is  reasonable  because 
the  transverse  stress  in  the  cracked  lamina  has  a  maximum  value  at  these  locations. 
The  process  continues  as  additional  strain  is  applied,  as  shown  in  Figure  32c. 

Reifsnider  et  ai,  in  Reference  [42],  suggest  that  when  a  matrix  crack  initiates,  the 
load  is  transferred  by  adjacent  plies  back  to  the  cracked  ply  over  a  distance  which 
depends  on  the  stiffness  of  the  cracked  and  adjacent  plies.  A  second  crack  will  form 
in  the  cracked  ply,  at  a  distance  away  from  the  initial  crack,  when  the  failure  stress  of 
the  matrix  is  exceeded.  This  process  continues  until  the  cracks  reach  an  equilibrium 
or  saturation  state,  commonly  called  the  Characteristic  Damage  State  (CDS).  The 
saturated  crack  density  is  independent  of  the  loading  history  but  depends  on  the  ply 
thicknesses,  orientations,  and  stacking  sequence,  as  discussed  in  Reference  [44]. 

A  number  of  studies  are  available  which  determine  the  stress  or  strain  for  the 
onset  of  cracking  and  the  spacing  between  cracks,  as  discussed  in  References  [43.  45]. 
.Some  studies,  such  as  References  [40,  41,  46],  report  the  changes  in  laminate  stiffness 
as  well.  Two  studies  on  laminate  stiffness  reductions  are  discussed  in  the  next  two 
subsections. 
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2.4.2  Stress-Strain  Behavior  of  Laminates  with  Matrix 
Cracking. 

Studies  on  the  stress-strain  behavior  of  laminates  with  matrix  cracks  were  performed 
by  Talreja  in  Reference  [40].  Talreja  determined  the  reduction  in  longitudinal  elas¬ 
tic  moduli  over  the  range  from  crack  initiation  to  crack  saturation  and  compared 
his  predictions  with  experimental  results  for  various  glass-epoxy  and  graphite-epoxy 
laminates. 

Consider  a  laminate  consisting  of  arbitrarily  oriented  plies  subjected  to  stress  in 
the  longitudinal  direction  (direction  of  0°  plies).  The  transverse  or  off-axis  plies  will 
develop  matrix  cracks.  Talreja  calls  these  plies  ‘cracking  plies'  and  the  remaining  plies 
are  called  ‘constraining  plies’.  Talreja  determined  that  the  initiation  and  development 
of  transverse  cracks  are  subject  to  constraints  from  two  sources:  the  thickness  of 
cracking  plies  and  the  stiffness  of  the  constraining  plies.  Talreja  reports  that  the 
strain  for  onset  of  matrix  cracking  decreases  with  increasing  transverse  ply  thickness 
and  increases  with  stiffness  of  the  constraining  plies. 

Talreja  uses  a  classification  scheme  to  qualitatively  assess  laminate  performance 
with  cracking  plies.  The  degree  of  constraint  is  classified  into  four  types,  as  discussed 
in  the  following  paragraphs.  The  laminate  stress-strain  behavior  corresponding  to 
these  four  types  of  constraints  is  shown  schematically  in  Figures  33a  to  33d. 

Type  A:  No  Constraint.  The  laminate  stress-strain  behavior  corresponding  to 
the  limiting  case  of  no  constraint  is  shown  in  Figure  33a.  The  response  is  linear  until 
the  longitudinal  strain  reaches  the  failure  strain  of  the  cracking  plies,  eppp*  where  ppp 
stands  for  first  ply  failure.  The  constraining  plies  are  unable  to  carry  the  additional 
load  shed  by  the  cracking  plies  and  extend  instantaneously  to  failure,  Cc- 

Type  B:  Low  Constraint.  The  behavior  of  a  low  constraint  laminate  is  shown 
in  Figure  33b.  Cracks  initiate  at  a  slightly  higher  strain  than  in  Type  A  laminates. 
Once  cracks  initiate,  the  cracks  multiply  under  constant  load  leading  to  an  abrupt 
increase  in  strain.  Further  increases  in  load  result  in  nonlinear  stress-strain  behavior 
due  to  the  formation  of  additional  cracks. 

Type  C:  High  Constraint.  Figure  33c  shows  the  stress-strain  behavior  for  a  high 
constraint  laminate.  Crack  formation  is  represented  by  a  knee  in  the  stress-strain 
curve.  Stable  crack  multiplication  requires  additional  load  increments.  Negligible 
nonlinearity  is  observed  due  to  the  low  increase  in  strain  with  matrix  crack  formation. 

Type  D:  Full  Constraint.  Figure  33d  shows  the  stress-strain  behavior  of  a  fully 
constrained  laminate  in  which  transverse  cracks  are  fully  suppressed.  The  laminate 
fails  when  the  failure  strain  of  the  constraining  plies  is  reached. 
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a.  Type  A:  No  Constraint. 

b.  Type  B:  Low  Constraint 
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Crack  Suppression 

c.  Type  C:  High  Constraint. 


d.  Type  D:  Full  Constraint. 


Figure  33.  Effect  of  constraining  ply  stiffness  on  stress-strain  behavior  of  laminates 
with  cracking  plies. 


Longitudinal  Strain  (percent) 

Figure  34.  Longitudinal  stress-strain  diagram  for  a  [0°,90|]s  glass-epoxy  laminate. 

Hahn  and  Tsai,  in  Reference  [39],  also  suggest  that  the  stress-strain  behavior  of 
laminates  depends  on  whether  the  test  is  load  controlled,  as  in  Figure  33b,  or  displace¬ 
ment  controlled.  In  a  displacement  controlled  test,  initiation  and  rapid  multiplication 
of  cracks  would  result  in  a  drop  in  laminate  stress  at  constant  strain.  This  type  of 
behavior  was  calculated  with  DYNA3D  and  previously  shown  in  Figure  14. 

Talreja  compared  the  measured  stress-strain  behavior  of  laminates  containing 
cracking  plies  with  the  behavior  predicted  from  a  damage  mechanics  theory^.  Con¬ 
sider  a  symmetric  cross-ply  laminate,  [0‘’,905]s.  made  of  glass-epoxy.  The  observed 
longitudinal  stress-strain  behavior  for  the  laminate  is  shown  in  Figure  34.  The  be¬ 
havior  corresponds  to  the  low  constraint  case.  Type  B.  An  abrupt  increase  in  strain 
and  subsequent  drop  in  longitudinal  modulus  occur  at  0.45%  strain.  The  reduction 
in  the  longitudinal  Young’s  modulus  with  applied  stress  is  shown  in  Figure  35,  where 
E\  and  Exo  are  Young’s  moduli  of  the  laminate  with  and  without  matrix  cracking, 
respectively.  The  observed  and  predicted  values  agree  quite  well.  Also  shown  in  Fig¬ 
ure  35  is  the  reduction  in  Young’s  modulus  predicted  with  a  ply  discount  scheme  in 
which  £22  and  G12  of  the  cracked  ply  are  set  equal  to  zero.  The  prediction  in  Figure  35 
based  on  the  ply  discount  scheme  overestimates  the  laminate  modulus  reduction  by 
12%. 

Additional  tests  and  analyses  reported  by  Talreja  indicate  that  the  laminate  stres.s- 
strain  behavior  and  the  reduction  in  laminate  modulus  depend  on  the  particular  layup. 

^The  predicted  stress-strain  behaviors  were  calculated  by  Talreja  using  a  damage  mechanics 
approach:  the  theory  characterizes  matrix  damage  by  a  vector  field  and  he  derives  relationships 
between  the  laminate  elastic  moduli  and  magnitudes  of  the  damage  vectors.  The  damage  vector  is 
orientated  normal  to  the  plane  of  the  crack.  The  magnitude  of  the  damage  vector  depends  on  the 
crack  density,  average  length  and  width  of  the  crack,  and  a  non-dimensional  factor  which  depends  on 
the  constraint  due  to  the  finite  thickness  of  the  cracking  plies.  Material  constants  for  the  undamaged 
state  and  one  damaged  state  (at  a  known  cra«;k  density)  must  be  measured.  With  his  analysis,  Talreja 
shows  that  an  initially  orthotropic  solid  becomes  anisotropic  as  a  result  of  matrix  cracking. 
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Figure  35.  Reduction  in  longitudinal  Young’s  modulus  with  applied  longitudinal 
stress  for  a  [0°,903]s  glass-epoxy  laminate. 


For  example,  the  behavior  of  the  cross-ply  laminate,  [0‘’,90°]s,  is  characteristic  of  the 
high  constraint  case.  Type  C.  This  is  because  the  alternating  longitudinal  plies  provide 
a  higher  constraint  to  transverse  cracking  than  those  in  the  [0°,903]s  laminate.  The 
ply  discount  prediction  overestimates  the  reduction  in  longitudinal  Young’s  modulus 
by  about  23%.  As  another  example,  Talreja  characterizes  a  [60°,  90®,  — 60°,  90°]2S 
laminate  as  approximately  a  no  constraint-type  laminate.  The  ply  discount  method, 
with  E22  =  0  and  Gu  =  0  in  the  90°  plies,  underestimates  the  longitudinal  modulus 
reduction  by  about  50%. 

Talreja  concluded  that  the  ply  discount  scheme  is  unreliable  for  predicting  min¬ 
imum  values  of  moduli  for  cracked  laminates.  He  further  concluded  that  the  ply 
discount  predictions  tend  to  overestimate  the  longitudinal  modulus  reduction  for 
high  constraint  cracking  laminates  and  underestimate  the  reduction  for  low  constraint 
cracking  laminates.  The  alternating  hoop  and  at-angle  ply  layups  typically  used  in 
motor  case  design  are  probably  high  constraint  type  laminates. 

Talreja’s  emphasis  was  on  reduction  in  the  longitudinal  modulus  of  cracked  lam¬ 
inates.  A  more  thorough  study  was  reported  by  Highsmith  and  Reifsnider  in  Refer¬ 
ence  [41]  in  which  reductions  in  Poisson’s  ratio  and  shear  modulus  were  also  reported. 
Their  results  are  briefly  summarized  in  the  next  subsection. 

2.4.3  Stiffness  Reduction  in  Composite  Laminates  with 
Matrix  Cracking. 

Highsmith  and  Reifsnider,  in  Reference  [41],  reported  results  of  quasi-static  tension 
tests  on  four  sets  of  glass-epoxy  laminates.  They  obtained  measurements  of  crack 
density  and  laminate  stiffness  verses  applied  longitudinal  stress.  Typical  results  for  a 
[0°,903]s  laminate  are  shown  in  Figure  36.  The  general  trend  is  that  as  crack  density- 
increases,  the  stiffness  decreases  until  both  the  crack  densifv  and  laminate  stiffness 
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Figure  36.  Crack  density  and  laminate  stiffness  reduction  versus  applied  longitudinal 
stress  for  a  [0°,903]s  laminate. 

reach  a  stable  value. 

Measured  reductions  in  laminate  elastic  constants,  reported  by  Highsmith  and 
Reifsnider,  are  compared  in  Table  6  with  those  calculated  by  APTEK  with  three  ply 
discount  schemes,  where  Ex  is  the  laminate  longitudinal  modulus,  uxy  is  the  major 
Poisson’s  ratio.  Gxy  is  the  laminate  in-plane  shear  modulus®. 

We  begin  discussion  of  Table  6  with  results  for  the  three  cross-ply  laminates.  The 
comparative  results  indicate  that  all  ply  discount  schemes  overestimate  the  reduction 
in  laminate  modulus,  Ex,  for  cross-ply  laminates.  The  largest  discrepancy  is  for  the 
[90§,0°]s  laminate  in  which  the  stiffness  reduction  is  overestimated  by  22%.  Signifi¬ 
cant  reductions  are  observed  in  Poisson’s  ratio;  all  ply  discount  methods  are  in  good 
agreement  with  the  observed  data  for  the  cross-ply  laminates.  The  measured  results 
also  indicate  that  that  there  is  little  (11.2%)  or  no  reduction  in  laminate  shear  modu¬ 
lus  for  the  cross-ply  laminates.  The  ply  discount  scheme  with  Gu  of  the  cracked  ply 
set  equal  to  zero  is  in  poor  agreement  with  the  observed  data.  This  suggests  that  it 
is  not  appropriate  to  discount  the  shear  stiffness  of  cracked  plies  loaded  in  transverse 
tension. 

For  the  laminate  with  45°  plies,  none  of  the  schemes  are  in  good  agreement  with 
the  measured  data  for  Ex,  or  Gxy-  Highsmith  and  Reifsnider  suggests  that  the 
agreement  is  not  good  because  other  damage  modes  were  observed  in  this  laminate, 
such  as  delaminations,  longitudinal  splitting,  and  fiber  breakage. 

2.4.4  Conclusions. 

Tests  and  analyses  results  presented  by  Talreja.  Highsmith  and  Reifsnider  indicate 
that  the  reduction  in  ply  properties  with  matrix  cracking  depends  on  the  ply  stacking 
sequence,  ply  stiffnesses,  and  cracking  ply  thickness.  Thus,  a  general  post-failure 

'’High.smith  and  Reifsnider,  in  Reference  [41].  also  reported  the  measured  and  predicted  changes 
in  laminate  bending  stiffness,  Dyy. 


Table  6.  Measured  and  calculated  reduction  in  elastic  constants  of  glass-epoxy  lami¬ 
nates.  Property  reductions  are  calculated  with  a  ply  discount  method. 

Laminate  Percent  Reduction  Discounted  Properties 

_ Ex  t/XY  Gxy _ 

[0,903)3  Observed  42.0  —  0.0 

Calculated  47.3  75.0  0.0  E22 

Calculated  47.3  74.3  0.0  E22  ^2\ 

Calculated  47.3  74.3  75.0  E22  ^2\  Gn 

[903,0]s  Observed  38^8  H^O  RO 

Calculated  47.3  75.0  0.0  E22 

Calculated  47.3  74.3  0.0  E22  i^u  ^2\ 

Calculated  47.3  74.3  75.0  E22  ^\2  ^21  G12 

[0,90]s  Observed  17.1  49.8  11.2 

Calculated  22.7  50.0  0.0  E22 

Calculated  22.7  48.9  0.0  F22  ^^12  ^21 

Calculated  22.7  48.9  50.0  E22  1^12  ^21  Gi2 _ 

[0,  ±45]s  Observed  9.4  8.1  28.4 

Calculated  3.9  20.1  9.7  E22 

Calculated  4.2  21.5  11.9  E22  I'u  ^2\ 

Calculated  30.1  -29.4  11.9  E22  v\2  ‘^21  C^i2 
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degradation  rule  must  be  implemented  into  DYNA3D  that  is  capable  of  modeling  ply 
stiffness  reductions  for  arbitrarily  stacked  laminates.  Allowing  the  user  to  specify  the 
degradation  behavior  for  each  mode  of  failure  would  make  the  code  more  flexible.  A 
flexible  post-failure  degradation  model  is  suggested  in  the  following  paragraphs. 

The  post-failure  degradation  rules  currently  implemented  in  DYNA3D  were  dis¬ 
cussed  in  Section  1.3.2.  For  the  fiber  breakage  mode,  all  the  lamina  elastic  constants 
{En,  E22,  1^12)  t/21,  and  Gu)  are  set  to  zero  over  an  arbitrary  number  of  time  steps 
(100  time  steps  in  the  current  version).  For  the  matrix  failure  modes,  some  of  the 
elastic  constants  are  set  to  zero  over  a  single  time  step. 

The  following  changes  to  the  current  DYNA3D  implementation  are  proposed: 

1.  Allow  the  user  to  specify  the  total  strain  increment  over  which  property  degra¬ 
dation  occurs  in  each  failure  mode.  The  simplest  model  would  be  a  linear 
degradation  in  elastic  constants  with  strain.  For  example,  if  matrix  cracking 
occurs  at  a  uniaxial  failure  strain  of  Ccrado  the  user  might  specify  a  reduction 
in  E22,  i^i2i  and  1/21  over  a  total  strain  increment  of  Ccrack/'f-  A  small  strain 
increment  would  produce  nearly  instantaneous  unloading  while  a  large  strain 
increment  would  correspond  to  constant  stress  after  ply  failure. 

2.  Allow  the  user  to  specify  the  percent  reduction  in  the  lamina  elastic  constants 
for  each  failure  mode.  For  example,  if  the  user  is  analyzing  a  cross-ply  laminate 
under  longitudinal  tension,  the  user  might  specify  an  80%  reduction  in  E221  t/u, 
and  i/21,  and  no  reduction  in  Gu- 

The  suggested  modifications  allow  the  user  to  model  constant  stress,  gradual  un¬ 
loading,  or  instantaneous  unloading,  following  tensile  or  compressive  failure  of  the 
fibers  or  matrix.  Material  property  reductions  as  a  function  of  strain  are  suggested 
because  laminate  stiffness  reductions  are  generally  reported  in  the  literature  as  a  func¬ 
tion  of  applied  stress  or  strain.  An  alternate  implementation,  to  the  first  modification 
listed  above,  is  to  allow  the  user  to  specify  the  number  of  time  steps  over  which  prop¬ 
erty  degradation  occurs  in  each  failure  mode.  The  simplest  model  would  be  a  linear 
degradation  in  elastic  constants  with  each  time  step,  as  is  currently  implemented  for 
the  fiber  breakage  criterion.  One  difficultly  with  reducing  the  moduli  over  a  specified 
number  of  time  steps  is  that  unloading  behavior  of  the  ply  will  not  necessarily  be 
linear  with  strain. 
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